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Abstract: We explore corrections to the "hard wall" gravity dual of QCD with two- form 
tensor fields bMN- These fields correspond to the quark bilinear = qa^'^'q of QCD, which 
generates states with quantum numbers 1 and l"* , or w/p-like and hi/bi-Wke mesons, 
respectively. We include new interaction terms, which render the model complete up to 
dimension six. We find that breaking chiral symmetry induces modifications to the spectrum, 
by mixing the vector current and , and breaking the degeneracy of the 1~~ and 1"' 
spectra. We also compute the hi — )• wvr decay rate and D / S ratio and the couplings of excited 
p states to two pions, performing a comprehensive new fit to the QCD data. We obtain a 
value for the magnetic susceptibility of the QCD quark condensate which is consistent with 
phenomenological estimates. These results allow us to evaluate the success of dimension 
counting and UV matching in this model. 
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1 Introduction 



Bottom- up AdS/QCD (pioneered in [1, 2] following earlier work on deconstruction [3]) suc- 
cessfully describes many features of low energy hadron physics. Bottom- up models do not 
descend directly from string theory (as do the top-down setups of [4, 5]), but rather posit 
the existence of a gravity dual for low-energy QCD, and study the dynamics of fields dual to 
QCD operators in a simple confining background. 

Though these frameworks are admittedly more ad hoc than top-down models, they remain 
essential tools for gaining a firmer grasp on the five-dimensional (5d) dual to QCD. For 
example, despite their greater predictive power and more "fundamental" origin, top-down 
models working in the supergravity limit represent a point in parameter space that is clearly 
distinct from the infrared (IR) limit of QCD.^ This supergravity (or a' — )• 0) limit relies 
on a large separation of scales between massless and massive string excitations. However, 
experimentally observed slopes of mesonic Regge trajectories suggest that no such separation 
exists: massive worldsheet excitations and massless ones all weigh in at ~ 1 GeV. Bottom- 
up models are not restricted to this a' = — >• limit, as they admit would-be massive 
string modes as new bulk fields.^ Bottom-up frameworks act as a simple testing ground for 
approximations made in top-down duals. Assuming that any bottom- up model we consider 
can be embedded into string theory, these simple frameworks may also shed light on the 
massive string sectors of top-down models. For example, as described in [7], the Lagrangian 
for the two-tensor field representing a leading massive string mode must be first order, a result 
which may shed light on the effective Lagrangian one would expect to derive from string field 
theory. 

Bottom-up frameworks have a significant disadvantage, however. Their Lagrangians are 
constructed term by term, instead of descending from a known string theory action: they 
contain a large number of undetermined parameters as a result. There is some disagreement 
in the literature about how to fix these free parameters (coupling constants, mass scales and 
operator scaling dimensions). For instance, [1, 8, 9] advocate fixing as many parameters as 
possible to the values they assume in perturbative QCD. The argument is that this procedure 
maximizes the predictive power of the model, leaving the fewest possible constants to be 
determined in the IR. 

We argue here that the above approach is somewhat counterproductive, as it might cause 
one to discard an otherwise viable framework simply because it undergoes renormalization 
group (RG) flow from the ultraviolet (UV) to the IR, and gives a poor matching to the low- 
energy spectrum evaluated using UV- valued parameters. It is important to remember that 
while the hard wall model lives on a space which is asymptotically AdS^, and thus has a 
conformal UV limit, the coupling constant is never small. The UV limit of the hard wall 
model is not the UV limit of QCD. One tacitly assumes that the UV boundary of the hard 

^Some attempts to go beyond this limit in top-down duals are explored in [6]. 

■^Relajcing the truncation at massless string modes and allowing the string length to be finite might also 
lead one to question whether the theory is local. In what follows we assume this to be the case. 
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wall model corresponds to an energy scale large compared to the scale of IR physics, but where 
the supergravity approximation of AdS/CFT is valid and the QCD coupling is already large. 
The hard wall framework thus approximates these intermediate energies with a "conformal 
window." As suggested by its original proponents [1, 2], though naive, this approximation 
works reasonably well for some quantities. Matching the UV behavior of this gravity dual to 
the UV of perturbative QCD presumes that correlation functions and scaling dimensions are 
not (or little-) affected by RG flow. There are certain operators - such as conserved currents 
- for which this should indeed hold true, but these represent only a fraction of the operators 
needed to describe the low-energy meson spectrum. We should treat all parameters not 
directly associated with conserved currents as free parameters in the bottom-up framework, 
and fit their values using low energy IR data. 

Note for instance that the extended hard wall framework suggested in [7] was subsequently 
analyzed in [9]. The authors determined all free parameters introduced with the two- form 
bMN by matching to perturbative calculations. The extension thus introduced no new free 
parameters. The holographic correlation functions of the vector and tensor operators were 
furthermore found to reproduce the large behavior expected from QCD. The resulting 
spectrum matched experimental data very poorly, however. Here, we leave as free parameters 
all quantities which are not explicitly protected under RG, and fit them to the low-energy 
data available. 

Even with a larger number of free parameters, AdS/QCD provides an appealing structure 
for modeling low-energy hadron physics that is motivated by top-down as well as phenomeno- 
logical reasons, and which naturally incorporates many previously-proposed organizing prin- 
ciples, like Vector Meson Dominance and Hidden Local Symmetry. Furthermore, by allowing 
these parameters to flow, one might even gain insight into the nature of the QCD renormal- 
ization group flow. 

In this work we continue our exploration of the hard wall model [1, 2] extended to include 
a two-tensor field dual to the QCD operator qa^^T°'q. With this extension, the hard wall 
framework includes all naive dimension three QCD operators. We focus our study on the 
vector (1 ) and axial vector (1"' ) modes generated by this operator, which should appear 
as massive worldsheet fields in a string description. We incorporate chiral symmetry breaking 
effects expected to appear at 0{Nc '^^'^). We also study the leading 0{N~^) modification to 
the mass of the 6i-like modes, performing a global fit to all free parameters in the model, 
including those of the original hard wall framework. For consistency, we add all interaction 
terms up to 5d dimension six. 

Our results shed light on the reliability of QCD dimension counting, the large A ap- 
proximation, and more. We find the interesting result that quantities associated with the 
massive string excitations tend to run significantly under RG, while most quantities related 
to massless excitations do not. 

The paper is organized as follows. In Section 2 we review the hard wall model and its 
extension to include all naive dimension three operators. We also describe the five-dimensional 
interaction terms that are responsible for chiral symmetry-breaking and leading corrections 
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of bi and p spectra. In Section 3 we review the requisite calculational techniques, relegating 
most details to the appendix. In Section 4 we present our results, including the deviation 
from perturbative QCD expectations. We conclude in Section 5. 

2 Model set-up 

We will now describe in detail the bottom-up hard wall model of [1] (and the extension pro- 
posed in [7]) that forms the basis of this work. Briefly summarized, the low-energy regime 
of QCD is characterized by the confinement of quarks and gluons to color singlet states, 
and by the breaking of the U{Nf)L x U{Nf)ii flavor symmetry of QCD to an approximate 
vector (isospin) subgroup. These elements are realized holographically via a gravitational 
background that induces confinement, and a non-trivial vacuum configuration of the quark 
mass operator qq. The other fields living in this background are dual to all dimension 3 quark 
bilinears, which generate mesonic states from the vacuum. On the gravity side the corre- 
sponding fields can be expanded in terms of an infinite (but discrete) tower of normalizable 
modes solving the equations of motion with different values of the 4d four-momentum (or 
meson mass). The vacuum state of the field dual to qq modifies the classical equations of 
motion for these fields, lifting the degeneracy between vector or pseudo-vector states. 

2.1 The extended hard wall model: background and fields 

The hard wall model lives on a slice of AdS^ truncated so as to induce confinement. We use 
the Poincare metric. 



where £ is the AdS^ radius. The rough equivalence between inverse radial coordinate z and 
energy scale allows one to identify with a UV cutoff which we eventually take to infinity 
(i.e. e — )■ 0). Meanwhile Zq^ corresponds to the scale of the effective theory we generate: the 
masses of the mesonic excitations are all proportional to Zq^. This somewhat artificial IR 
cutoff stands in for some more complicated geometry inducing confinement via a potential 
that goes gradually to infinity as 2; — )• (X) (see e.g. [10]). However, it is still an important toy 
model: in fact, some top-down duals, such as the Sakai-Sugimoto model, feature spacetimes 
that also end at a finite radial value. In addition to simplifying computations and making 
some of the IR characteristics of the theory more transparent, the hard wall model may thus 
shed light on some aspects of these top-down duals. Given this truncated geometry, we will 
have to make careful selections for the boundary conditions on the fields at the IR wall. We 
will discuss this in detail below. 

The fields living in this background are those dual to all of the naive dimension 3 quark 
bilinears of QCD as summarized in Table (1). These consist of the bifundamental scalar 
qq which corresponds to the complex scalar X, the left- and right-handed fiavor currents 
dual to left- and right-handed gauge fields of the gauge group U{Nf)L x U{Nf)ji, and finally 




e < z < zq 
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a two-tensor in the bifundamental corresponding to the complex two-form bMN- In these 
expressions a^y = ii'j^'ju — lulfj,)/'^, the adjoint index a = 1,2,3 for SU{2) flavor symmetry 
or a = 0, 1, 2, 3 for U{2), and the fundamental index i = 1, 2. 

operator dual AdS field 

g-4(l-7^)(?, iOf^=qiqj) Xi, (Re(X,,)) 
Ojl''' = q{x)t-^^lil+j')q{x) 
0^^'' = q{x)r^^\{l--i^)q{x) A-^j^ 
O^v = q{x)t''a^yq{x) blj^ 

Table 1. The operator/field correspondence in the extended hard wall model. 

It will prove most instructive and convenient to work with the conserved vector and 
broken axial-vector linear combinations of chiral currents, 

jV,a ^ oV,a ^ QL,a ^ QR,a ^ -^a^^^ ^ QA,a ^ QL,a _ QR,a ^ -^a^^^^^ ^ (2.2) 

whose dual fields are given by 

V=]^{Al + Ar), a = ]^{Al-Ar), (2.3) 

as these generate mass eigenstate resonances from the vacuum when chiral symmetry is bro- 
ken. For instance, the vector current generates the tower of J^^ = 1 or p-meson states, 
the divergence-free part of the axial current gives rise to a tower of 1^^ resonances, the ai 
states, and its divergence generates the = 0^ pion states. 

These operators and the scalar X have benefitted from extensive study in recent years. 
Our primary focus will therefore be on the tensor operator = q{x)a^vq{x), continuing the 
description of its spectrum and interactions initiated in [7]. As can be deduced from QCD, 
this operator encapsulates six physical degrees of freedom, and gives rise to two towers of 
resonances, one 1 and one 1"' . The latter are the bi resonances, nearly degenerate with 
the ai. As we will see below, these states also influence the existing hard wall spectrum as the 
"tensor /9s" mix with the "vector ps" generated by the vector current when chiral symmetry 
is broken. We should note at this point that this model is totally insensitive to isospin: the 
spectra of isospin- 1 and isospin-0 pairs (like the oj and p) are degenerate with each other. 
Because isospin-0 modes may undergo mixing with glueball states, which we neglect entirely 
here, we will focus our attention (except for a few special cases) on the isotriplet mesons. 

2.2 Action 

The action of our model supplements the original action for the hard wall model with a 
gauge-invariant first order action for the two- form field bj^iN'- 

S = S'hw + ^cs + Ssd + Sint , (2.4) 
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where 

5hw = |d^xV<7tr{|I)Xp-^|Xp-^(F| + Fl)} (2.5) 
5cs = ^ ^ (^5(^l) - cusC^k)) , (2.6) 

Here trF'^ = dto^. The additional interaction terms for hMN are described below. 

The action is fully invariant under the chiral gauge symmetry corresponding to the flavor 
symmetry of QCD; the bar denotes conjugation of the Lie algebra representation. The fields 
X, h transform in the bifundamental of U{Nf)L x U (Nf) The field strengths and covariant 
derivatives are 

{Fl,r)mn = dM{AL,R)N — Qn^Al^rim — i[{ALfi)M, {Al,r)n\ , 
DmX = dMX - i{Al)mX + iX{AR)M , 

DuhNP = dubNP - i{AL)MhNP + ibNp{AR)M ■ (2.8) 
We work with Hermitian generators and real structure constants normalized such that 

tr (f't^'^ = ^5-^^ , [e,t^] = if^H'^ . (2.9) 

Chiral symmetry breaking arises when we expand around a non-trivial vacuum of the 
complex scalar X, 

W - ^o(.) = ^ (2m,i---z--- + ^^f-^^") ^ ^ ' (^-^^^ 

which breaks U{Nf)L x U{Nf)R to the diagonal subgroup. Here A = 2 + \J + We 
choose the factors in front of the non-normalizable and normalizable modes in (2.10) such 
that ruq and a correspond to the quark mass and condensate, respectively. Since this result 
differs slightly from [1, 11], and since also the precise way in which the quark condensate 
a = {qq) is embedded in X will be important for us later, we briefiy review the derivation of 
(2.10) in Appendix A. The parameterization of Xq{z) in terms of the function v{z) is for later 
convenience. The original formulation of the hard wall model fixed the 5d mass of X such 
that A = 3, the bare scaling dimension of qq. This parameter should in fact run with RG, so 
we leave it as a free parameter of the model to be determined by comparison with low energy 
data. It is straightforward to see that evaluating X on its vev Higgs-es the axial-vector states, 
lifting their degeneracy with the vector combination of currents. 

The action for 6a/ tv above is written in terms of differential forms, with 

{Db)MNP = 3D[MbNP] ■ (2.11) 
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The fact that the action S^d is first order in derivatives is essential for constraining the 
number of physical degrees of freedom in bMN to match the six degrees of freedom contained 
in the dual QCD operator, while maintaining gauge invariance under the flavor symmetry. 
(See [7] for details). Naively, a massive complex two-form with a second-order action has 
twenty degrees of freedom in 5d. Eight of these, the fields, are Lagrange multipliers. 
If the equations of motion are first order, only six of the twelve remaining components are 
independent degrees of freedom. 

First order actions on spaces with boundaries require the addition of boundary terms to 
ensure a consistent variational principle [12]. This was first noted by [13] for free two-forms in 
AdS^. These boundary terms play a significant role, as without them the two-point function 
of the fields would vanish. In the hard wall model, there is also a boundary in the IR which 
requires a boundary term as well. However, it turns out that for our choice of boundary 
condition the on-shell value of this term vanishes. As there are no symmetries to guide us, 
choosing the appropriate IR boundary condition for Bmn may seem somewhat arbitrary, but 
in fact the choice is unique up to terms higher order in the fields, which are 1/Nc suppressed. 
In [7] we showed explicitly that of the two possible boundary conditions consistent with the 
variational principle, only one leads to physically reasonable behavior in the two-point 
function in the IR. 

Note the factor of sgn(/i) multiplying the bulk action in (2.7): as we allow fi to range over 
positive and negative values, this coefficient guarantees that the physical components of bMN 
have kinetic terms of the appropriate sign. While this factor is irrelevant to computations in 
the free 5d theory, the relative sign between the tensor kinetic term and the vector kinetic 
term has physical meaning once we couple these two sectors. 

In summary: the action described so far simply consists of gauge-invariant kinetic and 
mass terms for the fields in play. We have several free parameters: the scales i and zq, the 
couplings 55 and gt, and the AdS masses mx and /x. The only possible parameter not listed 
is a mass term for the fiavor fields, which we can set unequivocally to zero as it is related to 
the scaling dimension of the conserved vector current. 

2.3 Interaction terms 

Now let us turn to the interaction terms. Given that we work in a bottom- up framework, 
there are an infinite number of interactions we could include. We are most interested in 
exploring the effect of the new field bM n , particularly the effect of chiral symmetry-breaking 
in the quadratic order action (for which the most experimental data exists). In order to limit 
the number of new interaction terms, we include all terms only up to dimension six according 
to 5d dimension counting. Once the various fields are rescaled to have canonically normalized 
kinetic terms we find that they have mass dimension 

[X] = I [Lm] = [Rm] = I [bMN] = I . (2.12) 

The canonically normalized versions of the gauge fields have somewhat unconventional di- 
mensions due to the fact that the Yang-Mills coupling has dimension mass~^/^ in 5d. Using 



-7- 



this dimension counting, we find that there are a total of three interaction terms needed to 
complete the framework up to dimension six, the highest dimension present in the original 
hard-wall model. ^ 

First let us describe the terms which affect the quadratic order action for Bmn when 
chiral symmetry is broken. We begin by writing a term which produces the observed mixing 
between the vector meson states created by the vector current, and the vector mesons created 
by the tensor when chiral symmetry is broken. On the gravity side, this means that we 
should include a three-point interaction that involves the vector gauge field, buN, and the 
scalar X. When we expand around the non-trivial vev X = Xq, this contributes a quadratic 
coupling between Vm and Bmn- As discussed in detail in [7], there is only one such term 
allowed by the discrete symmetries of the theory with dimension less than six. This is a term 
of dimension 11/2, which we label with the coupling constant gi: 

Sg, =gij d^x^ii {bMNF^^X + XFl'^bMN + XF^'^Imn + Wi^'^^X} . (2.13) 

Chiral symmetry-breaking should also lift the degeneracy between the two towers of states 
(1 and 1^ ) created by the tensor . There is a unique dimension six term which serves 
this purpose. We label it with the coupling constant g2'- 

Sg,=g2i j d^x^ir {bMNXb^'^'X + bMNXW X] . (2.14) 

We have now described the interaction terms at or below dimension six, which contribute to 
the quadratic order action when chiral symmetry is broken. There is one additional term of 
dimension 11/2 not included in the gauge-covariant bMN action or the Yang-Mills action so 
far. We label it with g^ 

Sg, = igz I d^x^tr {6%F(^)p6^, + 6%F(^)p6^,} (2.15) 

and will find that it contributes significantly to the bi — t- lott coupling. 

Overall, we add three interaction terms, and three free parameters to the model: 

•S'int = Sg-^ + Sg2 + Sg^ . (2-16) 

Though we do not use Nc counting as our primary organizing principle for the interaction 
terms, especially from the perspective of the 4d theory, it is crucial to keep the Nc counting in 
mind. On the gravity side, working order-by-order in Nc essentially corresponds to working 
order-by-order in fields, as we now explain. Recall the standard large Nc counting rules for 
correlation functions of quark bilinears, which can include an infinite number of gauge fields. 

^More precisely, these are the highest dimension terms present in the 5d Yang-Mills action. Canonically 
normalizing the fields in the Chern-Simons action, however, one would find terms up to dimension 15/2. Most 
of these do not contribute to quantities of interest for us, however. 
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By counting factors of Nc that appear in planar correlators of r such masonic operators O we 
find 

(O1O2 • • • O.) ~ A^i"" . (2.17) 

We want bilinear operators to create mesons from the vacuum with amplitudes of order 1. 
We should then rescale these currents by factors of ^/IV^ such that 

((v^0i)(yiV^02) • • • iVN'cOr)) ~ iV]-"/' . (2.18) 

Reasoning analogously on the gravity side, we want to impose the requirement that the 
two-point functions of quark bilinears be order 1. This forces one to redefine the fields 
Al,m, Ar^m, bMN by absorbing appropriate factors of the couplings 55 and g^- With these 

— 1/2 

rescaled fields and the requirement that three-point functions should be order Nc , one 
can quickly deduce, for example, that 55 ~ 0{Nc ^'^^).^ Furthermore, after these rescalings 
we know the Nc order of any coefficient of an r-point interaction in the Lagrangian and can 
arrange terms in N~^ going order by order in fields. According to Nc counting with our 
rescaled fields, then, we have gigbQb ~ Nc gsg^g^ ~ Nc and g2g^ ~ ^c^- We allow 
these couplings to vary freely, fixing them by fitting to IR data, and checking whether their Nc 
scaling approximately matches expectations. Note that modifying any of the terms we have 
written above by adding derivative operators would certainly increase their 5d dimension, but 
would leave them at the same order in Nc'- it is truncating at dimension six that allows us to 
remain with a finite number of terms. 

2.4 Gauge fixing 

Some brief remarks about gauge fixing are warranted, as we choose a different gauge from 
[1]. The pion and sigma fields, with the pions being the Nambu-Goldstone bosons of chiral 
symmetry breaking, are naturally encoded in the fiuctuations in phase and magnitude of X, 
respectively. While the vacuum breaks the U{Nf)L x U{Nf)ji gauge symmetry, the theory 
itself still possesses the full symmetry. The simplest way to determine the physical spectrum 
is to make gauge choices that fully fix the gauge symmetry. Note that we work with the 
5d theory on shell and in the leading saddle point approximation: we are thus dealing with 
classical gauge fixing. 

In our model we only consider pion fiuctuations, and not sigma modes. Thus, X takes 
the restricted form 

X(x,z) = Xoe^^-^^'") , (2.19) 

where it{x, z) = tt^(x, z)T"', with tt" real. We can use part of the gauge symmetry to eliminate 
the pion fiuctuations from X. Since X transforms as a bifundamental, the action of the gauge 
group is 

X ^X' = e-'^^Xe'^^ . (2.20) 

* While clearly the order 1 coefficients of these terms should run with RG, we assume that their leading 
order in Nc remains unchanged between UV and IR. 



- 9 - 



By taking = — = vr, we have that 



X' = Xo . (2.21) 

This is our first gauge fixing condition. We will drop the prime on X from now on. 

Since X oc 1, the gauge transformations leaving X invariant are the diagonal U{Nf)v 
transformations generated by A^, = = Ay. Under these the axial vector field A transforms 
as a charged matter field, while the vector gauge field V transforms (not surprisingly!) as a 
gauge field. We will use the remaining gauge freedom to impose conditions of V, while A is 
allowed to be a general 5-vector, which we parameterize as 

AM = iAj^ + d^ip,A,) . (2.22) 

Here d^Aj;^ = 0, and we may assume d^d^ip 0, as it would otherwise have been included in 
Aj^. The modes of are dual to the tower of axial vector mesons, while (p, A^ contain the 
modes dual to the pion tower. 

A gauge transformation can always be found that brings us to = gauge. The 
remaining gauge freedom is to make U{Nf)v transformations that preserve this condition; 
these correspond to z-independent gauge transformations. Writing = + d^4'v, where 
d'^V^ = and d'^(j)v 7^ 0, one finds that the equation of motion for (py following from (2.4) is 
d'^dz4>v + 0{(^'^) = 0, where $ represents any field. Working perturbatively, this means that 
(j)y has a leading order z-independent piece plus corrections that are quadratic in the fields. 
We may use our final bit of residual gauge symmetry to set the leading order piece to zero. 
It can be shown that the quadratic order terms in (py do not contribute to any three-point 
couplings when plugged back into (2.4). Hence for the purposes of this paper we may safely 
take 

Vm = {V^.,Q), d^'V^ = Q. (2.23) 
3 Computational techniques 

We wish to understand the effect of the chiral symmetry breaking interaction terms on low 
energy observables such as masses of the mesonic excitations, their decay constants, and 
some interactions that mediate observed decays. In this section we describe in general terms 
how to extract this information from the extended hard wall model. Details are left to the 
appendices. 

We begin by analyzing the quadratic part of (2.4), which we denote S^'^\ The linearized 
equations of motion and boundary conditions that result from imposing = win deter- 

mine the bulk-to-boundary propagators and spectrum of normalizable modes. From these we 
derive the two-point functions, which contain the meson masses and decay constants. Hav- 
ing properly identified the physical fluctuations, we then compute a selection of three-point 
couplings. We apply some of these results to obtain the tree-level contribution to the decay 
rate and D/ S wave amplitude ratio for hi — )• wk. 
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3.1 Equations of motion and boundary conditions 

We would like to fix X to its vev (2.10), and study the action (2.4) on this background as a 
functional of the 5d fields V, A, b. One should immediately raise objections to this program 
since our slice of AdS space together with the scalar vev do not constitute a solution to the 
coupled Einstein-matter equations of motion. However, provided that the mass of the scalar 
is in a certain range, there will exist a solution to the coupled equations with asymptoti- 
cally AdS boundary conditions [14], (see also [15]), and our treatment serves as a reasonable 
approximation to that exact solution. The mass range is —4 < < 0, corresponding to 
2 < A < 4. The lower limit on the mass-squared is the Breitenlohner-Freedman bound, 
below which an asymptotically AdS solution becomes perturbatively unstable^. The upper 
limit is easier to understand from the dual field theory point of view. A > 4 corresponds 
to an irrelevant operator. Deforming the field theory by an irrelevant operator renders it ill- 
defined in the UV; the corresponding gravity statement is that asymptotically AdS boundary 
conditions can no longer be maintained when > 0. This conclusion is also backed up 
by a perturbative analysis of the coupled Einstein-matter equations along the lines of [16]. 
Henceforth we will restrict ourselves to the range 2 < A < 4 and work in the approximation 
where backreaction of the scalar vev is neglected. 

Variation of the action at quadratic order results in the linearized equations of motion 

n -WA , ^ ^ A '^9i9b sgn(Ai) . . 929^ sgn(/x) . . 2 ^ r 1 a 

U = idb+ — -kb "^-^ v(z)-kdV viz) -kb, (3.1 

Q = di.{dV -^v{z){b + b)^ , (3.2) 

{) = di.dA-4,v{zfi^A. (3.3) 

The variation reduces to a set of boundary terms when evaluated on a solution to these 
equations, namely, 

55(2) ^ j ^4^^^ I l^^l^^.upa (^5T ^ _ ) + ^(a.^^ - ^^,A,)5A^'+ 



+ 



2 ^ , 2giv{z) T ^ ^° 



z 



(3.4) 



Here is the self-dual (-I-) or anti-self-dual (— ) part of 6^1/, which can be obtained via 
projection, 

^% = \ [s^^K^ ± \^,u'") &P<x . (3.5) 

The top sign in (3.4) is chosen when /i > and the bottom sign when /i < 0. See [7] for 
details. 

^We could extend the range of A down to the unitarity bound A > 1 by choosing A to be the smaUer of 
the two roots of A(A — 4) = m\, which is permissible when —4 < m\ < ~3. 
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Boundary conditions should be chosen such that the UV boundary terms at z = e and 
the IR ones at z = zq in (3.4) vanish: this guarantees that we are in fact expanding around 
a stationary point of the action functional. We discuss each of these in turn, beginning with 
the UV boundary terms, for which we will require some results concerning the asymptotic 
behavior of solutions to (3.1)-(3.3). 

The equations of motion for V, A are second order, b, meanwhile, has two first order 
equations involving both the real and imaginary parts, which we can convert into a single 
second order equation by acting with the appropriate differential operator. As usual, the 
general solution for each field will be a sum of two types of solution with differing behavior 
as 2; = e — )• conventionally referred to as normalizable and non-normalizable. The leading 
scaling behaviors are 

non-normalizable: b^i, ~ z^'^' , V^,^^ ~ -z" , 9^ ~ z'' , 

normalizable: b^^ ~ z^"^ , V^,Aj^r^z\ ~ ^mm{2,2A-4} _ (3 

These results follow from a detailed analysis of (3.1)-(3.3) given below. Some comments and 
caveats should be stated here, however: 

• When // > (< 0) it is (bj^jy) that contributes the dominant behavior in (3.6). 

• The scaling (3.6) is independent of gi,g2- The interaction terms, being proportional to 
the tachyon vev, give subleading terms in the differential equations as z — )• 0. Therefore 
these results can also be obtained in the 91,52 — ^ limit where we have the original 
hard wall model and a decoupled free 6-sector. Note that these statements are only true 
provided A < 4 which we are assuming. 

• We assumed a further restriction on parameter space when stating (3.6). The gi terms 
in (3.1), (3.2) couple V with specific components of b. This means that each solution 
for V induces one for b and vice-versa. There are thus a total of four solutions to 
the coupled b-V system, with correlated asymptotics b ~ {z^^^^ , z^~^ , z^~^) and V ~ 
(^^t\iJ'\+g-A ^ z*^, z^). In order that the V wavefunction induced from the non-normalizable 
fe-type solution not dominate the wavefunction of the non-normalizable V-type solution, 
one must assume |/i| < 6 — A. The necessity of this condition can be argued from both 
physical and technical standpoints. First, the leading asymptotics of V are directly 
tied to the dimension of dual vector current, qj^f^q. This is a conserved quantity; its 
dimension is protected and given by the classical value Ay = 3. On the gravity side this 
requires the leading asymptotics V ^ z^. Secondly, as we show in Appendix C, allowing 

> 6 — A leads to divergences in the on-shell action that cannot be removed by 
counterterms that are local functionals of the sources; the holographic renormalization 
procedure of [16, 17] breaks down. 

The two classes of UV boundary conditions we consider are differentiated by whether or 
not we turn on non-normalizable solutions. If we are only interested in normalizable modes of 
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the 5d fields, dual to meson states, we set the coefficients of the non-normalizable solutions to 
zero. This eliminates the UV boundary term in (3.4) in the e — )• limit. If we are interested 
in the bulk-to-boundary propagators of the 5d fields, with which we can study correlation 
functions of QCD operators, we turn on non-normalizable solutions. The UV boundary term 
in (3.4) is made to vanish by holding the fields fixed, 5$ = 0, at the boundary. The values 
they are fixed to are the classical sources dual to the QCD operators. The precise boundary 
conditions, given in 4d momentum space, are 

Reb%ik,e) = ^^i^r^,(A:) , V^ik,e) = Vl^{k) , 

Aj^{k,e) = Al{k) , ip(k,e) = 7:(k) , (3.7) 

where ^{k, z) = f (i^xe~*'^'^<I>(x, z), and the self-dual or anti-self-dual projection on b is again 
correlated with the sign of All other components of 5d fields are determined in terms of 
these asymptotics via the equations of motion. The sources 7^jy(/i;), V^{k), and A^{k)+ik^7r{k) 
are dual to the (classical) dimension-three tensor, vector, and axial-vector quark bilinears, 
respectively. In particular, n^k) is dual to the divergence of the axial current, which contains 
the same information as the scalar quark bilinear. 

Next, let us consider the IR contribution to (3.4). The boundary term is a sum of terms 
of the form P$ • 5$, where P$ is some function of the fields. Naively, the two options are 
a Dirichlet condition that holds $ fixed, or a Neumann condition that sets P$ = 0. In the 
case of 6V^ and 6A^^ we can use U{Nf)L x {7(A^j)/j gauge invariance to argue in favor of the 
Neumann condition, as the Dirichlet condition would explicitly break this symmetry. Since 
bfj,u does not correspond to a conserved current in QCD, there is no analogous argument to 
determine its boundary condition. In [7], however, we used physical considerations to argue in 
favor of the Neumann condition. In the free limit, gi,g2 — ^ 0, choosing the Dirichlet condition 
results in a pole at fc^ — >• in the tensor-tensor two-point function. As there are no massless 
particles in QCD with the appropriate quantum numbers to explain this divergence, it is 
reasonable to conclude that this boundary condition does not yield a good gravity dual to 
QCD. With the Neumann condition there is no such divergence. We continue to follow this 
prescription here, since we will be interested in regions of parameter space where 51,52 are 
small. This results in the following set of IR boundary conditions:^ 



9195, 



= , 

z=zo 



b^^ik, zo) = 0, d,V^ + '-^v{z){b^, + 6^,) 

dzAjil=,, = , {d,^ - A,)l=,, = . (3.8) 

These conditions apply whether we are restricting to normalizable modes or determining the 
bulk-to-boundary propagators. 



^As we discussed in [7], one can modify the IR boundary conditions by adding terms to the action (2.4) 
that are localized on the IR boundary. Such terms, beyond the ones already considered in (2.4), are subleading 
in 1/Nc, and are on the same footing as additional 1/A^c-suppressed terms that we could add to the bulk. 
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While the 5d-covariant notation of (3.1)-(3.3) provides a relatively simple presentation in 
which symmetries are manifest, the dynamics are obscured by the redundancy of the descrip- 
tion. For example, 6a/ iv contains twenty real components, but only six of these correspond 
to physical degrees of freedom. Since we will only be using on-shell quantities in the 5d 
theory, we are free to use equations of motion to eliminate extraneous fields. This process is 
summarized here; see appendix B for further details. 

A minimal set of fields may be chosen as follows. First, Vf^,Aj^ contain all degrees of 
freedom associated with the vector- like uj/p and axial /i/oi mesons respectively. Below we 
will often use rescaled versions of these fields, 

Vf,{k, z) = g^V^{k, z) , Aj;{k, z) = g^A^,{k, z) . (3.9) 

The scalars ((/?, Az) contain all t^/tt fluctuations. The linearized equation of motion for A^ has 
no z-derivatives; thus Az can be solved for algebraically, leaving a second-order equation in 
dz for V9(/c, z). We will find it more convenient, however, to introduce a different field variable 
in the pion sector. We define A{k^ z) such that 

^{K z) = - , .^ dzA{k, z) , Az{k, z) = -^^A{k, z) . (3.10) 

KZqU(z)^ Z^V[Z)^ 

The relation between ip, Az implied by these identifications is consistent with their equations 
of motion. Here k = yfW is always real and positive since we we will only be using these 
expressions on shell. 

Of the twenty real components of huN^ eight are eliminated by a Proca-like condition. 
Half of the remaining twelve have the interpretation of momenta and half coordinates. Using 
the equations of motion one can show that 6a/ at may be parameterized in terms of two 4d- 
transverse vectors, H^^ik^z) and YJ^ik^z), as follows: 

h,AK z)= - '-^{e,rk, {n. + ^m + ^^^^dzk[,n,^+ 

fcZo [ |/U| +U2 ' 

2sgn(^)z 2{kzo)ui \ 

-Ozki^W^] + -r-, K[^V^] '> , (3.11) 



IH - U2 ' ' |/i| - U2 

b,z{k,z)= -g,V-A T\'^'\ w, + ^OzV,- ^7,^f > 4 . (3.12) 

[Zo{\n\-U2) \fJ'\-U2 2:0(1^1+1*2) J 

Observe that these expressions also have dependence on V^. This comes from mixing terms 
in auxiliary equations of motion that have been used. We have introduced notation for 
combinations of parameters that appear often in the following: 

ui{z) = 2gigbg5v{z) , U2{z) = g2glv{zf . (3.13) 

As is clear from (3.11), H parameterizes the real-transverse components of b^u, while W pa- 
rameterizes the imaginary-transverse components. The transformation properties of 6 under 
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charge conjugation and parity identify H with the mesons and W with the tensor- hke 

w/p's. 

The factors of (|/i| it 7x2)"^ arise from ehminating auxiUary components (such as h^z) 
and plugging back into the equations of motion. As we will discuss soon, we must restrict 
parameter ranges such that > |ti2(2;)|, Vz G [0, zq], otherwise we will not necessarily obtain 
a positive definite spectrum. Intuitively, the 52 couplings act like off-diagonal terms in the 
mass-squared matrix for h once X is evaluated on its vev. If these off-diagonal terms become 
too large, the eigenvalues may become negative. Since v{z) is a monotonically increasing 
function, the constraint is equivalent to > |ti2(2;o)|- Note that this restriction also makes 
sense from the perspective of the large expansion. The combination of couplings 525^ that 
appears in U2 should be 0{N~^) compared to the diagonal terms of the mass-squared matrix, 
which, by construction, are 0{1). 

In appendix B we show that the quadratic part of (2.4) is equivalent, on shell, to the 
following action, written in terms of fields {A, Af_i,'H^, Ri_i = (W^, V^)"^}: 

S^^^ = - I -0^ dzti |u;(^)l(A;2 - oW)i + w^a)Mk' " 0^'^^)^'^+ 

+ ^^(fe)^M(^' - O^'^)^^ + i?MW(^)(A;2i _ 0'^P))R^^^ + Sb„dry , (3.14) 
where all operators are of Sturm-Liouville type. Explicitly, 

O^""^"^ = - [^Md.] + , (3.15) 

OC-) = - wT^.d. + - ul) , (3.16) 



z 



The weight metric in the p sector is 

W(p) = diag(w(^),'u;(^)) , (3.18) 

and the weight functions are 

z I f u\ \ 

^ ' z^d/il - ^2) ^ z V - U2J 

Here we see that positivity of the norm in the 61 and p sectors imposes the conditions \p\ > 
\u2{zo)\ and |/i| > U2{zo) + Ui{zo) on the parameters of the model. The last term in (3.14) 
encompasses all of the boundary terms present and determines the on-shell action to quadratic 
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order in sources when we evaluate S^"^^ on the bulk-to-boundary propagators. We will have 
more to say about this below. 

Finally, let us state the boundary conditions satisfied by these fields. First we introduce 
some notation for this purpose that will come in handy below. In each sector, we define linear 
functional, Q and P, which encode the UV and IR boundary conditions, respectively: 

Qin)[f] = -4^in)dJ , Qia)[f]=f, [/] = 3 (/ " ^O^W^^/) ' 

Q{p) = {Q(w),Q{v)f , with 



and 



P{n) [/] = /, P{a) [/] = ziw^a)dj , P[b) [/]=/ + ^0^(6)5^/ , 

P{p) = {P{,w),P{v)f , with 

P{w)[if,9) ] = / + ZoVJ{w)dzf 7-J-2 9 , 

l/^l ^2 

P{v)[{f,9) ] = ZoW^y)dzg 7-. / . (3.21) 

IMI ~ ^2 



Additionally we define sources 



1 .n... „ 



S{p)ij. = iS(yj)iJ,^ , with 

Letting {$} = {A,An,'Hn,Rfj,} denote the collection of field content for all of the different 
sectors and O the corresponding Sturm-Liouville operators, the bulk-to-boundary propagators 
are solutions to {k^ — 0)$ = 0, satisfying the UV and IR boundary conditions 

Q[mk,e) = S, P[^]{k,zo) = . (3.23) 

Here Q, P and S refer to the appropriate boundary conditions or sources for a given field 
in {$} . Normalizable modes corresponding to individual meson states are eigenfunctions of 
the various O's. In this case we set Q[^]{k,£) = 0; as e ^ this isolates the normalizable 
modes. The IR boundary conditions P[$](A;, zo) = thus lead to quantization of the spectra. 

3.2 Bulk-to-boundary propagators, on-shell action, and two-point functions 

Although wc must ultimately resort to numerics to obtain explicit solutions for bulk-to- 
boundary propagators and normalizable eigenfunctions, a great deal of progress can be made 
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analytically in evaluating the on-shell action and two-point functions. The on-shell action at 
quadratic order is given by the boundary terms in (3.14). Our analysis in appendix B results 



m 



'S'bndry 



(2^ 



tr 



+ 



1 



+ 



w 



{v) 



sgn{fj,)kui 



sgn{iJ,)kzoUi 



|/i| - U2 



+ 



^tr {pM(A;,e)-Q[cI>](A:,e)} 



(3.24) 



Here we have used the IR boundary conditions to eliminate the boundary terms at z = zq. 

We would like to evaluate this action on the bulk-to-boundary propagators. The Q's can 
simply be replaced with the sources using (3.23), but in order to express the P's in terms of 
sources we need to find the explicit solutions for the bulk-to-boundary propagators. Since we 
will eventually be taking the e — )• limit, it is useful to separate out the non-normalizable 
and normalizable components of the solution. 

3.2.1 TT, ai, and hi sectors 

We begin with the vr, oi, and bi sectors. The non-normalizable and normalizable solutions, 
f^-'^-{k,z) and f^\k,z), are two particular solutions to (/c^ — 0)f = distinguished by 
their asymptotic behavior as z — >• 0. The general solution is a linear combination / = 
gii.n. jn.n. _|_ ffi. jn. ^ The Coefficients are fixed by the boundary conditions (3.23) and this leads 
to the bulk-to-boundary propagators for vr, ai, and bi as 



<i>{k,z) 
C{k) = 



r-^-{k,z)-C{k)r-{k,z) 



Q^-^ik,e) 
^"■°-(fc,zo) 
P^ik,zo) 



C{k)Q^ik,e) 



S{k) 



where 



(3.25) 



Here we are using a shorthand Q'^ "' = etc. Acting on an element of {$} with the 

corresponding Q and evaluating at z = e leads to the UV boundary condition, while acting 
on it with P and evaluating aX z = zq gives zero, thanks to the form of C{k). Plugging this 
solution into (3.24) gives a contribution 



Pmk,e).Qmk,e) 



P--{k,e)-C{k)P-{k,e) 
^^'''Q---ik,e)-C{k)Q-ik,ef^''^ ' 



(3.26) 



3.2.2 p sector 

Essentially the same formula holds for the p sector, suitably generalized. In the p sector 
(fc^ — 0)R = is a coupled second order equation which has a four- dimensional solution 
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space. We isolate four particular solutions by their behavior near z = 0. They are non- 
normalizable and normalizable solutions of w- and w-type. A solution of "tu-type" can be 
characterized as one that becomes pure W, R = (W,0)"^, in the (71 — )• limit where tensor- 
vector mixing is turned off. Thus we may also denote these as tensor-like and vector-like. If 
i G {n.n.{w),n.n.{v),n.(w),n.(v)} is an index running over these four solutions, then a general 
solution has the form R = c*i2*. 

The coefficients are fixed with the four boundary conditions Q[R] = S, P[R] = 0. Let a, (3 
be indices running over {w,v}. Define the two-component vectors c^''^', c^' with components 
(.n-'^-iP) gQ on. In the p sector the boundary conditions are now given by two- by- two 
matrices Q"^-"-, Q"^-, P'^-"-, P'^- with components {Q''-''-)afi = (3(a) [-R''-"-^^^] etc. Then the 
boundary conditions take the form 

Q"-'^- {k, e)c'^-°- + Q^- {k, e)d^- = S , P"-'^- {k, zo)c°-'^- + P"' {k, zoK = , (3.27) 

the solution to which is 

d^- = -C(A;)c°-^- , c°-"- = [Q"-'^-(A:,e) - Q'^-(/t, e)C(/t)]-^ 5 , with 

C{k) = {P-ik,zo))-'p----{k,zo) . (3.28) 

This gives the coefficients c*, and hence the bulk-to-boundary propagator c*i?* in terms of the 
source. In order to compute the on-shell action we are interested in -P[-R], which takes the 
form 

P[R] {k, e) = P"-^- {k, e)c^-°- + P'^- {k, e)c"^- . (3.29) 
Plugging in the solution for the c's, we find that the p sector contribution to (3.24) is 

m(t, e) ■ Qimt, e) = 5W . (3-30) 

using the p sector versions of P, Q and S, as usual. Here T denotes transpose and the matrix 
in the denominator is to be understood as matrix multiplication by the inverse from the left. 

3.2.3 On-shell action and holographic renormalization 

The contributions (3.26), (3.30) determine the on-shell action at finite e. In order to compute 
field theory correlators via gauge/gravity duality, we need to take e — )■ 0. Generally there 
are both divergent and finite terms in the limit. The finite terms can be further divided 
into those which are entire functions of k'^ and those which are non-analytic in k"^. In all 
cases considered here the divergent terms come from the S{Q^''^')^^ P'^'^' S piece of (3.26), 
(3.30), and are furthermore entire functions of k"^. Finite terms may come from this piece as 
well, but they are also analytic functions of k'^. The non-analytic finite terms are obtained 
from the S{Q^''^')~^C{k)P'^' S piece. The Q"' term in the denominator of (3.26), (3.30) never 
contributes in the e — )• limit. These results are reviewed in Appendix C, where we present 
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the necessary asymptotics and, in particular, give our conventions for the non-normahzable 
and normahzable solutions. 

The process of holographic renormalization [16, 17] introduces a counterterm action to 
cancel the divergences that arise in the on-shell action. This mechanism is dual to the 
regularization of UV divergences in the field theory. The counterterm action should be a 
local functional of the sources, as guaranteed by the analyticity in momentum space of the 
terms that diverge as e — >• in the on-shell action. Similarly, the finite terms in the on-shell 
action that are entire functions of k"^ are scheme-dependent; they may be shifted to zero by 
introducing local finite counterterms. The renormalized action is thus determined by terms 
finite in e and non- analytic in k. The results of Appendix C lead to 

to, r (fk r mlk'^e^-^CMik) 2C(a)(k) , 



/2|Mh2 

+ - — 



2\^Ji\ 



per 



+ ^'"^'"^^if + C,^m k^V^'vA , (3.31) 



where Cap are the components of Cp. This expression is valid for generic parameter values. 
In the special case A = 3, the non- normahzable solution in the pion sector has logarithmic 
instead of power law behavior for small z and the prefactor in the first term should be replaced 
with (6 - 2A) 1. 

3.2.4 Two-point functions 

This result can be used to compute two-point functions of the QCD quark bilinear operators. 
On the one hand, the sources appear in the generator of QCD correlation functions as 

W=l^exv{i j d^x (yOA*«0^'« + (A^^" + 8^)0^'" + T^'"'Oj^i'')^ ^ , (3.32) 

where the expectation value is over the QCD path integral. On the other hand, gauge/gravity 
duality tells us that W = Zgj-g^-y[V^ , , ir ,T], the gravity partition function with asymptotic 
boundary conditions given by the sources. 

This relationship must be taken with a grain of salt, of course, as we are certainly not 
working with a well-defined theory of quantum gravity. We could imagine that our model has 
an embedding into string theory, and the action (2.4), evaluated on shell, gives the leading 
contribution to Zgrav in a gs and a' expansion, Zg^av = e*"^ + • • • . These correspond to the 
1/Nc and low energy expansions in field theory language. However, the string dual of QCD is 
not known, and in any event the usefulness of these expansions (in particular the low energy 
expansion) for QCD is questionable. Until these issues are improved upon, the approach of 
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AdS/QCD should be understood from a different point of view: we use Z^s = e*'^ to define the 
generator of correlators in an effective theory, W^s = Z^s , and we see how well this effective 
theory describes QCD. 

The on-shell action (3.31) is sufficient for determining the two-point functions in our 
effective theory of QCD. From the relationship Wes = e^^ we schematically have (O1O2) = 
— (5^i(502e*'^|<^(,=o, where c/iq denotes the collection of sources. The results are conveniently 
expressed in a set of matrix elements 11 which are gauge invariant Lorentz scalars, defined in 
terms of the (momentum space) correlators as follows: 



{0^'''{k)0^^\-q)) = i(27r)^5(^)(fe - q) [{k\, - k^K)U^^^{k) + k^k^U^^^k)] 6^' , 
{0^''^{k)0^'\-q)) = i{27r)'S^^\k - q){k\, - A;^fc,)n^'^(A;)<5'^^ , 

(Oj^(fc)0,?(-g)) = ^{27:)'6(^Kk - q) {v^.,,.Il^'''Hk) +rl,,.Il^'^Hk)) 6^" , 
{0'j:fik)0];'\-q)) = (27r)4<5(4)(fc - q) (k.VuX - k,r^,x) n^'^(fc)5'^'' . (3.33) 

The transverse and longitudinal tensor projectors are defined as 

i'p^r.i = k%6:^ - 2k[^k^X] ' = ^k^^k^x] > (3-34) 

and satisfy = k^V. Making use of these definitions, some straightforward computation 
yields the following set of two-point functions: 

f2\^l\-2 /2|/x|-2 



11^'^ {k) = 2 1,2 2^'"'"^k) ) Yl^'^ iji) = \7a+i i^wvik) + Cvw(k)) . 

95'^ ^0 2gbg^kz^Q 

(3.35) 



All of the dynamics is contained in the functions C{k), which must be obtained numerically. 
They are given by C = (pn.-j-ipn.n. ^ evaluated at the IR boundary. The precise definitions 
of the normalizable and non-normalizable modes for each sector can be found in Appendix 
C. 

3.2.5 An interesting one-point function: magnetic susceptibility of the quark 
condensate 

We can also use (3.31) to compute one-point functions of operators in the presence of external 
fields, yielding, for instance, the magnetic susceptibility. 

Let us briefly review how the magnetic susceptibility of the quark condensate is defined 
[18]. Consider turning on an external magnetic field — a classical source whose electromagnetic 
potential we denote A^^. In such a background one expects the tensor current Oj^lf" = qf^a^iyq 
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to acquire a vev, and the magnetic susceptibility measures the ratio of the tensor current vev 
to the external field. More precisely, we have a susceptibility Xi for each species of quark qi, 
defined through 

{qicrf,uqi) = eiXi{mi)Ff^ , (no sum on i) , (3.36) 

where Cj is the electric charge of the quark, and F^J^ = dfj^A'^ — d,yA'^^. To be concrete we will 
consider the case of the up quark, = 2e/3: 



{ua,,u) = {Oil' + Oli') = |x«(0^'° + 0''')F^i . (3.37) 




The holographic prescription together with the on-shell action (3.31) can be used to 
compute X- The electromagnetic current of the quarks 



can be expressed in terms of fiavor currents J^'"" = qt^-j^q according to J™ = e(| 7^'*^ + J^'^). 
Thus by choosing our vector current sources V^'" to be of the form 

V^'' = , V^'^ = eA^^ , V^'^ = V^'' = , (3.39) 

we find that the exponent of (3.32) contains the term VJ?'"gt"7'^g = (A'^^)^J^'^. This shows 
that A'^^ is an electromagnetic background field. Now we evaluate the tensor one-point func- 
tions appearing in (3.37) in the background (3.39) with all other sources turned off. To leading 

(2) 

order in the electromagnetic field we can use W = exp (iSren), and functionally differentiate 
with respect to the tensor source. We find 



^ _sgn(/x) 

r=Ao=TT=o 



. , 2 iCu,vik) + C,^{k)) {fv)% , (3.40) 



where % = 5'fe(^o/^)'''' ^ and (/y)" = -2i/cr f"'", leading to 



{Oli\k) + Oli^k)) = lllf^l'^l {CUk) + C^k)) F^i . (3.41) 



In principle we could study a spacetime-dependent magnetic susceptibility, but in order to 
keep things simple and compare with the literature later, we will take the A; — )• limit to 
recover the constant background value. 

Notice here that we are working perturbatively in the background field. To calculate the 
one-point function in a finite background field one should expand the 5d fields around the 
constant finite background solution = V^, (3.39), much as we expand around the nonzero 
tachyon vev, (2.10). This type of calculation was recently carried out in [19] in the context 
of the model proposed in [7], but with parameters fit to the UV as in [9]. Computing the 
tensor and scalar one-point functions in this background, the authors of [19] find a x(^'^') to 
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all orders in the external field. We are computing the linear response, x = x(0)- We will 
compare with their results as well as other theoretical and lattice results in section 4.2. 

Since we are computing the linear response to the field it is sufficient to evaluate the 
scalar vev's on the right-hand side of (3.37) at their leading-order ^^'-independent value. 
This result is given simply by the quark vev: 

(O^'O) = ^(Q-g) = I and (O^'3)=0. (3.42) 

(See Appendix A for details). Combing the k ^ Q limit of (3.41) with (3.42) and plugging 
into (3.37) we find 

4sgn(^) 1 
ly ■ hm - 

In section (4.2) we give numerical estimates for this result based on our best-fit parameter 
values. 



Xu = „ I _ 2 • 1™ Z (Cn,v{k) + a^(fe)) . (3.43) 



3.3 Normalizable modes, masses, and decay constants 

The physical quantities we ultimately wish to extract from solving the linearized equations 
of motion and from the on-shell action are of course the masses and decay constants of the 
resonances. There are two equivalent procedures for deriving this information: (1) from the 
pole structure and residues of the two-point function in the form given in (3.35), or (2) by 
writing the two-point functions explicitly as a sum over momentum space poles whose residues 
are given in terms of the tower of normalizable modes describing the bulk fields (as described 
in more detail below). It turns out that the first form makes most apparent such general 
properties as the /c^ — t- oo or /c^ — )■ behavior of the correlator, which the infinite sum form 
conceals. The second repackaging, meanwhile, becomes more useful when we are interested 
in specific (usually light) resonances and their decay constants, which, in turn, the first form 
leaves obscure. 

In order to derive the sum-over-poles-type expression, we should first discuss the towers 
of normalizable modes. In each sector define the L^-normalizable Hilbert space of functions, 

'H = {(l)eLl[e,zo]\QmKe) = Pmk,zo) = 0} ■ (3.44) 

One can show that each O is an Hermitian operator on the Hilbert space Ti of a given sector. 
In the vr, ai, and bi sectors O is furthermore positive definite, while we expect O in the p 
sector to be positive definite for fixed momentum and gi sufficiently small, though we do not 
prove it. 

Given the form of (3.14) it is clear that the eigenvalues of the various O's should be 
interpreted as the mass-squared of normalizable modes. The O's of the 7r,ai, and bi sectors 
are independent of the momentum k"^, so we get a standard dispersion relation for the 4d 
states, A;^ = m^, where m = 0, 1, 2, . . . labels the eigenvalues. In the p sector we see that O 
depends on k. It can be shown that the eigenvalues depend only on the momentum-squared, 
and thus we get a nontrivial dispersion relation k'^ = Am(fc^)- The p masses are the solutions 
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i/v7(MeV) 




to this equation; we denote them by k'^ = m^. In Figure 1 we plot the momentum dependence 
of the first three eigenvalues \m- The p masses are given by the intersection points of the 
curves with the diagonal. Notice that when the momentum gets large enough the eigenvalues 
become negative; the momentum-dependent off-diagonal terms in O, (3.17), are competing 
with the diagonal terms. We work with parameter values such that the momentum scale at 
which this occurs is well beyond the scale where we expect our description to be valid. Near 
each mass eigenvalue there is an expansion 



Xm{k ) = m„ + 



2 dXm 



d(A;2) 



ml) + O {{k'^ - m'if) . (3.45) 



We have an orthogonal basis of eigenfunctions for the Hilbert space of each sector, T-L, 
satisfying 

0(pm{z) = m^cpmiz) for the iT,ai,bi sectors 

04)m{k, z) = Xm{k'^)(t>m{k, z) for the p scctor (3.46) 



and normalized such that 



dzw(j)m4'i = ^mi for the vr, oi, bi sectors 



d2;(/)^W(p)0; = 6ml for the p sector. (3.47) 

We have chosen our basis functions to be real-valued for convenience. Note that (j)m in 
the p sector is a two-component column vector of functions. When we want to specify its 
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components we write 0" where a, /3 run over w, v as usual. The bulk Green functions for the 
operators (A;^ — O) can be expressed in terms of the normalized eigenfunctions as 

G{k, z; z') = ^rniz)4>i{ ) ^ sectors 
^-^ — m£, 

m '"- 

r. n '^ ST <Pmik, z)(p^ (k, z') 

Gai3{k, z;z) = 2_^ — p — — sector. (3.48) 

They satisfy, in each sector, 

(/c^ — 0)G = w~^6{z — z ) for the vr, oi, 61 sectors 

{k^ - 0)G = W-^(5(2 - z') for the p sector. (3.49) 

We can use these bulk Green functions to obtain convenient expressions for the bulk-to- 
boundary propagators by generalizing the argument in [20]. First observe that the following 
Green's identities hold, as can be explicitly verified using (3.15)~(3.17), (3.20), (3.21). Let 
V'(z),x('2) be arbitrary functions on [e,2:o]- Then 



dzw [xiy - O)^ - - 0)x] 



for vr, ai, h\. (3.50) 



For the p sector let be the components of a two- vector and matrix- valued 

function. Then we have the Green's identity 



£ dz [xo. (W(A:2 - 0))„^ V'T/? - ^^7/3 (W(A;2 - O))^^ x.] 



-^0"' {QoXxWomn - {QomfiPo.{x\) 



(3.51) 



Here the notation {P(^p-^ct)\^\)p means that P acts on the first or second column of the matrix 
'(/', with the result indexed by /?. Now employ these identities as follows. Take x = ^ to be the 
bulk-to-boundary propagator, and "0 = G to be the Green function. On the left-hand sides the 
second terms vanish because the bulk-to-boundary propagators are solutions to the equations 
of motion for arbitrary fc^, while the first terms integrate to give $(^')) using (3.49). On the 
right-hand sides both terms vanish at the IR boundary because both satisfy Neumann 
conditions there, while only the •P\G\ terms survive on the UV boundary since Q\G\ = 
there. Noting that Q[^] is the source 5, we find the desired relations: 

$(/c, z) = 4>m{z) for the vr, oi, 61 sectors. 



'^(^'^) = E^^Sl2^^4M'^™(^'^) for the sector. (3.52) 
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With these expressions the boundary action, (3.24), takes the simple form 

{7r,ai,6i} V ^ rn. ™- 



1 f d^k ^ ^\P[(t>m]{k,e)-S{k)\^ 



(3.53) 



p sector 



This is the sum-over-modes form of the on-shell action, vahd at finite e. Our next task is to 
determine the renormahzed action in the e — t- limit. In particular, we want to determine 
how the -P[0] for each sector are related to the normalizable modes, f^ ., which we can compute 
numerically. 

Let's begin with the vr, ai, and bi sectors. In each of these sectors, the (j)m are eigenfunc- 
tions of O and thus solve the same equation as the bulk-to-boundary propagators, but with 
k"^ —7- . The general solution is a linear combination of normalizable and non-normalizable 
modes, f^-{mm,z) and /'^■^■{mm, z). Thus UV boundary condition Q[(j)m] = fixes the 
relative normalization such that 

0(m^, z) = N (rimm, z) - ^^^^I^r^-^m^, z)] . (3.54) 

Imposing the IR condition determines the eigenvalues, while N in each case is fixed by the 
normalization condition (3.47). As e — )• 0, one finds from the asymtotics in Appendix C 
that the second term vanishes relative to the first for any finite z. Hence, in the limit, the 
eigenvalues are given by 

P[/'^'](mm, zq) = P'^'inim, zq) = for mass eigenvalues nim ■ (3.55) 

This indeed corresponds to the poles of the C{k), (3.25), that parameterize the two-point 
functions (3.35). Meanwhile the normalization constants are given by 

N-^= wrimm,zf + 0{e>'') . (3.56) 
Jo 

Using the explicit aysmptotics of the /'s we can also compute 

PMi&ie) = iV(.,,r.) (1 + 0{e>^)) , P^a0Z>]{e) = 2%,^) (l + 0{e>^)) , 



P(,)[C^](e) = 2Ar(,^^)^i^ (1 + O(5>0)) . (3.57) 

^0 

The p sector requires more care due to the momentum dependence of the eigenvalues. 
Eigenfunctions are solutions to the equation (A — O)0 = 0. This equation is slightly different 
than the equation satisfied by the bulk-to-boundary propagator since A must be treated 
independently of the momentum k'^ appearing in the p sector operator O. Nevertheless it is 
easy to show that the leading asymptotics of the four linearly independent solutions are the 
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same so we will continue to use the same notation introduced around (3.27). The general 

solution is R = ^^^c^R', and in this case the UV boundary conditions Q[R] = may be 
written as 

Q'^-"-(A,/t,e)c'^-°- + Q"-(A,A;,e)c°- = . (3.58) 

Using the asymptotics of the Q matrices in Appendix C, we find that the c^'^' are linear 
combinations of the c^' with coefficients that go as positive powers of e. Hence for finite z 
we can neglect the contribution of the non-normalizable modes to the eigenfunctions 0^^^ . In 
particular the IR boundary conditions which determine the eigenvalues become 

P"-(A,A;,zo)c''- = detP°-(A,A;,zo) = A G {A™(A;2)} . (3.59) 

If we go on shell by setting Am(^^) = k'^ then (3.59) corresponds to the poles of C(p)(A;), 
(3.28). 

Notice that for each Am in this sector, the ffist equation of (3.59) determines the ratio 
of the normalizable mode coefficients appearing in Let us denote these coefficients = 
w,m) 1 N{v,m))'^ SO that we have 

P''-{Xmik^),k,zo)-Nmik) = . (3.60) 

This equation does not determine the A^'s outright, since the rank of P is one. The overall 
scale is fixed by the normalization condition which, as e — ?• 0, becomes 



1= / dz{^mfW^p)^„ 

Jo 





(3.61) 

where in this sector (Am(fc^), ^, -z) = {f(uim)'9^wm)^'^^ ^^'^ terms that vanish as 

e — )• have been neglected. The A^'s, like the eigenvalues, depend on the momentum. In 
order to determine the on-shell action in the p sector we also need 

P[4>m]{k, e) ■ S{k) = S{kf ■ P- {Xm{k^), k, e) ■ 

^ ^13^it^)^N^^^^)k-%,{k) + -N^,^^)V^{k) , (3.62) 

as e — )■ 0. 

We can now collect results and evaluate the on-shell action, (3.53). We find that the 
leading terms as e — )• are finite; naively there are no divergences. One may wonder how this 
is consistent with the analysis in section 3.2, where renormalization was necessary. However 
in this section we have been analyzing the limit mode by mode, tacitly assuming that it 
commutes with the sum over modes. This is clearly false, taking into consideration (3.52): this 
expresses something that finite in the limit, the bulk-to-boundary propagator, as an infinite 
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sum of terms that vanish in the hmit. In the following we assume that commuting the limit and 
the sum is equivalent to passing to the renormalized action, at least for computing the physical 
quantities of interest — masses and decay constants. We have checked numerically that the 
formalism of the previous section and this one give the same results for these quantities. Thus 
we find 



lim Si^l 



e-5>0 



1 



^0 



tr < A; vrvr 



E 



Nf , 

\iT,m) 



— m? ^ 



+ 



95 



Nf , 

(a,m) 



+ 



H > T' 



■fiv 



Ar2 

{w,m) 



V. 



Ar2 



A;2 



m 



(6l,m) 



+ 



Nf , 

(v,m) 



+ 



8i sgn(/i)£l'^l-i 



jpa^ 



{k')S 



(3.63) 



where we have introduced the additional notation Vj^ = k'^rj^'^ — k^k^ . 

Varying with respect to the sources, applying the holographic prescription, and taking 
into account the matrix element definitions (3.33) leads to 



^4 



Nf , 

('!T,m) 



" m 



n 



A,A 



m 



/i^2 ^4 



E 



(a,m) 



A;2 



(ai ,m) 



n 



T,T± 



2^2|m| 



-2 



„2r.2 2|m|+2 



E 



A;2 



(bl,m) 



n 



5(|A;2z,^ 



E 



m '^^(p.m) 



n 



n 



TT\\ 



2im 



-2 



A;2 



2|m|+2 
-0 



E 



(w,m) 



{p,m) 



(fc2) ' 



2sgn(^)^l^|-l ^ 



(3.64) 



These are essentially the sum-over-poles expressions for the matrix elements that we can 
compare to field theory templates to read off decay constants. See for example [7]. (To 
compute residues for the 11 involving p spectra one still needs to use the expansion (3.45). 
Also the momentum-dependent Nj^p-^ should be evaluated at the pole.) We find 



f{m.) 

J TT 



eVm 



pTm 



N 



(7r,m) 

2 ' 
''^{7r,m)^0 



rm 
J ai 



2N 



(a,m) 



h 



bi 



2(^/^0 



2A'"(i,,m.)("i(p,m)) 



V^N^w,m){m(^p,m)) (£/~o)l^l~l 

9bm(^P,^)Z^^J'^- >^[p,m)i^t, 



(3.65) 



,m) > 



where the prime indicates differentiation with respect to /c^. Let us summarize. Expressions 
(3.65), together with (3.55) and (3.59), constitute our final analytical results for masses and 
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decay constants. We have boiled them down to quantities that involve simple manipulations 
with the normalizable-mode wavefunctions, f^'. Masses are obtained by imposing the gener- 
alized Neumann condition at the IR boundary, P[f'^']{zo) = 0. Decay constants are expressed 
in terms of the normalization constants, N = \\f^'\\~^, where the norm is the natural one on 
the Hilbert space T-L, (3.44), (for e — )■ 0). The normalizable-mode wavefunctions in turn are 
easily computed numerically. 

3.4 Three-point couplings 

Having exhausted the quadratic-order Lagrangian, we now predict three-point couplings that 
are either new to the hard wall model (such as the one mediating bi — )• coir), or modified due 
to the mixing between the Imn sector and the vectors {gp-Kn)- For the purposes of comparing 
to experimental data, we are only interested in tree-level couplings of the lowest-lying states. 

States in field theory are dual to normalizable modes in the bulk. We associate to each 
state \ sn) a 4d field Sn{k), and we identify these fields with Fourier coefficients of corresponding 
bulk fields expanded in the eigenfunction basis: (p^^^k^z) = J2n^n{k)(pn\z). The holographic 
4d effective action for the fields Sn{k) is obtained by plugging these mode expansions into the 
bulk action and KK-reducing along z. The mode expansion for each bulk field is summarized 
in Appendix D. By design, these expansions diagonalize the quadratic part of our model: 



K^k^ - <n))^"'^ - - A(p,n)(fc'))p"'^ . (3.66) 



Due to the nontrivial dispersion relation in the p sector, the canonically normalized fields are 

pl{k)=\l-\\.{ml.)X'^~pl{k) . (3.67) 



We now turn to the computation of some three-point couplings. 
3.4.1 5b. 



Let us begin with the vertex that generates bi — )• wvr (or equivalently hi — )• pir). bi modes 
are contained in 6a/ iv and vr modes in Am-, while p modes can be either in 6a/ at or Vm- Thus 
we are looking for terms in the bulk action that contain either two factors of 6a/ Af and one 
Am, or terms that contain one each of 6a/ at, Am, and Vm- In principle, S^a, Sg^, and Sg^ 
can all contribute, but it turns out that Sg^ only contains three point couplings of the form 
b-V-V or b-A-A. Hence the entire coupling comes from the gauge covariant derivative in the 
bMN kinetic term, or from the g^, interaction term. The relevant contribution from the kinetic 
term is 

_isgi^ j tr{6AD6-6Am} 
^ " ^4^ / {{hMNhpQ + hMNbpQ)An] , (3.68) 
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while the relevant portion of the interaction term is 

Sg, D -253 J d'x^tr { Im [6^m6%(9^Ap - dpA^")] } . (3.69) 

We must take Ar — )• Afft"- with e SU{Nf) C U{Nf) to guarantee that we get a pion mode. 
Meanwhile, one of the b factors should live in SU{Nf) (to give the bi) and one in U{1) (to 
give the w). 

Inserting the mode expansions (see Appendix D), we eventually obtain after a straight- 
forward but tedious calculation the following three-point couplings among excitations in the 
b, uj, and vr (with momenta k, p, and g, respectively): 



_ 2zo 



bmnrP ' Q ~l" CmnrQ 

•fc)6'^''"(fc)-w"(p)+ 



TT'^^^iq) . (3.70) 



Hered^Q = {2TT)-^S'-'^\k+p+q)d'^k-d'^p-d^q is the measure obtained from Fourier transforming 
a product of three fields to momentum space, and we have used the notation uj for the U{1) 
component of the p field, (3.67). The dimensionless couplings a, b, c, are given by 

i^OimnribifYinri Cmnr) — ; ' {dmnribmnriC-mnr) i (^•'^-'-) 

with 

"i(7r,r) r° 



J'mnr 



_ sgn(//)m(^,„)Zom ^(^^^^^^^| 



+ 



|/i| - U2 

~ m ~ 

bmnr 



'm(uj,n)'m(^^r)ZO Jo [ 



+ 



+ 93gbghw^b,)W(^^)&^zd,4>t^ (4-) + ^^''^zod,<t>^:A I . (3.72) 

Note that since our framework is insensitive to isospin, Shp-w is identical to the above with 
6" — )• /i, a; — )• p". The first subscript in each of the above coefficients refers to the mth state 
in the bi tower, the second to the nth state in the uj tower, and the third to the rth state 
in the vr tower. We are primarily interested in the first (lightest) state in each tower, but in 
theory this framework allows us to compute all of them. 
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With this 4d action in hand, we now compute the decay rate 61 — )• wvr, as weh as the 
D/S ratio, which characterizes the angular dependence of the amphtude. We consider the 
lowest mode in each tower, defining a = am, and similarly for 6, c. 

3.4.2 hi decay rate 

The matrix element for the decay process, in terms of a basis of orthogonal, transverse 
polarizations for the massive vector-like modes, becomes 



with 



A = 



2zQ 



/ 



{—ak ■ p + bp ■ q — cq ■ k) , B 



2zq 2 



/ 



m^{a-b-c). (3.73) 



The indices i,j = 0, zt, label all possible polarizations for the 61 and for the uj. From here we 
find the 61 decay rate via standard techniques. The absolute square, averaged over incoming 
61 polarizations and isospins and summed over outgoing u polarizations and tt isospins, is 
found to be 



1 



A' 2 + 



jp ■ kf 
ml ml 



{k-qf {p-qf , {k-qfip-q) 



ml^ml mlml 



+ 



■mlmlm^ 



+ 



2AB (1 - ' - ' + 1^ ■ ^^^^ ■ ^)(^ '-P) 



ml^ml mlml 



ml mlml 



Working in the lab frame, we consider momenta of the form 

k^, = {mb^ , 0) , Pf, = iV'ml + ct^, -q) , = (\/m2 +q2, q) 

The decay rate is 



bl—^OJTT 



|qp^ ( mfei - + |q|2 - y^m^ + |q|2) 



Snmb^ Jq \/m2^T]qFv'm|+l^ 

01 

where is the root of the argument of the (5-function, 

1 



4ml^ 



( I I 2, 
[rrib, -m^-m^ 



2\2 



2 2 



m\ 



bi 



(3.74) 



(3.75) 



(3.76) 



(3.77) 



The decay rate for the hi is identical except that we don't need to average over initial state 
flavors; thus rhi_,p^ = (iV| - l)Tbi-^^-;,. 
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3.4.3 D/S ratio 

In addition to the polarization-averaged decay rate, we can predict the angular profile of 
the final state. This information is encapsulated in the D/S ratio, or the amplitude ratio 
for D-wave versus S-wave decays. The ratio is defined from the amplitudes via some basic 
quantum mechanics. The initial 61 state has spin 1, and no angular momentum. The final 
state must thus have J = 1. Since the tt has 5 = and the u has 5 = 1, they may have some 
relative angular momentum in the final state, either L = 2 or L = 0. Parity forbids L = 1. 
The relation between the initial and final states is given by the Clebsch-Gordan coefficients. 
For angular momentum L = 2, we have 

|1, l)j ^ yf >^2,2|1, -1) - y^l2,l|l,0) + /^^2,0|1, 1) , 

|i,o)j ^ y^y2,i|i,-i) - y|i2,o|i,o) + y^y2,-i|i,i> , 

|1, -1) J ^ Vf ^2,-2|l, 1) - /^^s.^ill' 0) + \/^^2,o|l, -1) , (3.78) 

where the orbital angular momentum of the final state has been written in terms of spherical 
harmonics and the spin has been written as \s,ms). If L = 0, then the 5 quantum numbers 
must be the same as the J's: 

|l,m)j^yo,o|l,m) , (m = 1,0,-1). (3.79) 

The general final state produced in 61 decay can be written as a linear combination of these 
possibilities that depends on the initial spin of the 61 : 

\i,i)b, = 5yo,o|i,i).. + ^ 

11,0)6, = 5yo,o|l,0)^ + D 

\i,-i)b, = 5yo,o|i,-i).; + i) 

By squaring the above eigenstates, we can see that the coefficients 5, D should satisfy 
|5p + |Dp = 1, but their ratio is determined by the dynamics of the theory via the matrix 
element (3.73). 

We define the polarization-dependent amplitude A4^^ generically as 

M'' = {l,j\u,T\l,i)b, . (3.81) 

So 

T|l,l),, = A^++|1,1)^+7W+0|1,0)^ + 7W+-|1,-1)^ , 
T|1,0);„ = A^o+|l,l)^+A^o°|l,0)^+A^°-|l,-l)^ , 
T|l,-l)fe, = M-+\l,l)^ + M-^\1,0)^ + M—\l,-l)u^ . (3.82) 



y2,2|l,-l) 





/fl2,o|l,0) 




^2,-2|l,l) 




y2,o|i,-i) 



;3.80) 
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Evaluating the amplitude (3.73) for various choices of polarization vectors allows us to deter- 
mine D and S. For this purpose we need an explicit basis of polarization vectors. We work 
in the rest frame of the bi, so the polarization vectors corresponding to spins ms{bi) = 0, ±1, 
are 

' (^) ^ ~ V2 ±i 

We can boost an identical set of rest-frame polarizations to arrive at the polarization of u 
which moves in the lab frame with three-momentum — q: 



(3.83) 



-(0) , 



1 



\ 



1 



I 



1 



\ m^{m^+E^)1z1± J 



(3.84) 



where qx,Qy,Qy are the components of the three- vector q, and q± = qx i: iqy. Note that 
e ■ q = 0, with given in (3.75). Also we have = y^mj + |qP. 

As an example, consider the evaluation of Ai^^. Taking (3.73), and plugging in the 
explicit four-momenta and polarizations, leads to 
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±± 
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1 fE„ 



2\m. 



1 1 sin^ I 
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ITlh-^ 2 ■ 2 , 

2-q sm I 



a-1{a(v^-i) + b^.,}(i- 



cos 



+ 



1 /I67r 
6 



|yl(7r+^-l) +S^xy|y2,o • (3.85) 



In the first line we introduced standard spherical coordinates in q-space, and in the second 
line we defined 

x=\q\/muj, y=\q\/mn- (3.86) 
Analogous computations yield the following set of matrix elements: 
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(3.87) 
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where S and D are given by 



S 
D 



3 



A + + +B 



xy 



A 



1 + x2 - 1 



B xy 



(3.88) 



Note that we should not expect these D,S to satisfy |5p + |Dp = 1, as we computed the 
matrix element only to leading order in perturbation theory. At this order, our prediction for 
D/S is 

Af^/^T^-l) +5^xy 



-V2- 



+ x^ + 2]+B'^xy 



(3.89) 



3.4.4 Qp^^ 

The mixing of the hMN field with the vector field implies that not only the mass spectrum, but 
also the interactions of the 4d p meson are modified. We must therefore re-compute the p-n-jr 
coupling which appears in [1], and include the revised version when we fit the undetermined 
parameters of the extended model. 

Since the pions are contained entirely in Am while the /?'s are contained in both Vm and 
bMN, we are interested in bulk three-point couplings of the form V-A-A and b-A-A. It turns 
out that only S'hw and Sg^ contribute to this coupling, with 



"S'hw ^ 



^91 ^ 



abc 



(fx 



dz 
z 



X I ^v{z){-{d''K^)ip'd'^ip^ + 2hl,Ald^^'^} , 



(3.90) 



where Hmn = Re 6a/ at- 

We can now sum these two contributions, go to momentum space, and plug in the mode 
expansions of Appendix D. The result should be compared with the standard definition. 



Sp^n = gpn.r'^ I d''xpl{x)Tl\x)d^'Tl'{x) 



(3.91) 



We find the following on-shell value for the coupling: 



9. 



95 



pTTTT 



^/l^\r^im^~) Jo 



■ _ sgn(^)m(,,^)ni 

zo(|^|-n2) 



(3.92) 



One can verify that this reduces to the result of [1] in the gi,g2 ^ ^ limit. 
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4 Results 



Based on the above computations and those of the original hard wall model, we have made 
predictions for the following physical observables: nip, nib, nia, rUn, f^, /J, fb, fa, fn, Qp-Kir, 
and the D / S ratio and decay rate for hi — )• ujtt. In principle we have KK towers associated to 
each of these, but to compare with data we restrict nip and gp-j^^^ to the first three p states, m^r 
to the first two pion states, and everything else to the lowest KK mode. These observables 
depend on the free parameters of the original hard wall model, zq, 55, m^, a, and A, and 
on the new parameters we introduced by including the AdS two-tensor field hMN and its 
interactions: gb, ffi, 52, 53 -"^ 

We allow all parameters to vary, instead of fixing them to values calculated in perturbative 
QCD. This stands in contrast to the original hard model [1] and to [9] where the interactions 
proposed in [7] were also included: these works fix as many observables as possible to values 
computed in perturbative QCD. Of course, since the vector flavor symmetry is preserved, the 
vector current acquires no anomalous dimension. This implies that we can fix the 5d masses 
of the corresponding gauge fields to zero. All other parameters do vary between the UV and 
IR, however. It is more accurate (though admittedly less predictive) to fit all parameters by 
fitting to IR data. As we have a large number of free parameters, our main purpose will be 
to explore how well AdS/QCD models conform to hadronic data when we push them beyond 
a' = and Nc — )• 00. 

As described in detail above, the tensor field 6a/ 7V and its chiral-symmetry-breaking 
interaction terms induce mixing between the vector mesons generated by the QCD two-tensor 
O^, and the vector mesons generated by the QCD vector current. These terms also lift the 
degeneracy between the vector- and pseudo- vector states generated hy huN- In other words, 
coupling this new sector to the original hard wall model modifies predictions for observables 
like nip and gp-K-n', furthermore, the predictions we make for the previously neglected bi mesons 
depend critically on the parameters of the original model. In order to make consistent and 
accurate predictions, therefore, we perform a comprehensive fit to the experimental data 
including all free parameters. 

The full parameter space is subject to some restrictions, discussed in subsection 3.1, and 
summarized here. 

• We restrict 2 < A < 4, where A is the scaling dimension of qq dual to the scalar X. 
Values of A less than 2 violate the Breitenlohner-Freedman bound, while it was shown 
in [16] that A > 4 for the classical solution of X ruins the asymptotically AdS behavior 
of the geometry. 

• We restrict < 6 — A. This guarantees that the V^-type solution induced by turning 
on a source for bpi, does not dominate the Vp mode induced by turning on a source for 
Vp itself. Violating this restriction implies UV scaling behavior for Vp at odds with the 

'^The AdS radius (. does not play a role for these observables but will for the magnetic susceptibility. See 
section 4.2. 
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scaling dimension of the conserved vector current Jy, which is protected from RG flow. 

> 6 — A also leads to divergences in the on-shell action that cannot be removed by 
local counterterms. 

• We must also impose some constraints guaranteeing that the Sturm-Liouville problems 
defining the masses and decay constants of each sector remain well-defined. For certain 
values of the free parameters, the weight functions appearing in the Sturm-Liouville 
inner products pass through a pole and switch sign. To guarantee that this does not 
happen, we require constraints 



Finally we note that the observables we have computed only depend on the relative sign of 
ji and gigbdn- The fit prefers them to have the same sign, but it does not determine whether 
this sign is positive or negative. 

4.1 The global fit 

The experimental and lattice data we use for the fit are summarized in Table 2. We have 
chosen the data which we expect to be most accurately modeled by our results. For instance, 
isospin and 1 states are degenerate in our framework, so we have chosen to fit our predicted 
masses to the measured masses of isospin 1 states (i.e. the p, the bi and the oi), as the isospin 
modes may mix with glueballs. Furthermore, the error bars in the D/S ratio are given by 
the average PDG estimate, but these errors may in fact be much larger when systematic 
effects are included. For example, the D/S ratio is found, experimentally, to have a non- 
trivial phase - something we cannot reproduce in our framework. It is suggested in [25] that 
the phase may be due in part to final state interactions between the vr and the oj, which we do 
not take into account. Similarly, the gp-j^ couplings generically assume complex values [26], 
though in our mode they are real by definition. However, the experimental data indicates 
that the phase is close to either 1 or —1 (within the error bars), so we fit to a single positive 
(or negative) number, instead of considering phase and magnitude separately. 

Values for the parameters which participate in the fit are given in Table 3. The overall 
rms value of the best fit result is 0.296 it 0.026. The quantities marked with square brackets 
in Table 2 are not included in the global fit. fb-^ and fj are the only two observables which 
depend on gb alone (not in combination with gi or (72), while r(6i — >• cjtt) and the D/S 
ratio are the only two quantities which depend on g^. The values quoted in the table are 
found by using fj to fix gb, which we then use to determine fbi- These quantities are clearly 
incompatible in the existing model.® Meanwhile, we fix g^ using the value of r(6i — t- utt) and 
use it to compute a prediction for the D/S ratio. 

In order to facilitate comparison to the hard wall model, (Model B of reference [1]), we 
define combinations of certain parameters that are dimensionless in both frameworks. The 

*If we were to perform the reverse operation, fixing (?(, using the /hi we would find /J ~ 13.7 MeV. 




and 
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quantity 


fit 


experimental (lattice) result 


reference 


p 


777.1 ± 5.0 


775.49 ± 0.34 


[211 

L J 


rrio' 


1433 ± 6 


1465 ± 25 


[21] 


P 


1783 ± 11 


1720 ± 20 


[21] 




141.7 ±0.8 


153 ±7 


[22] 


[fl 


184* 


(184 ±15) 


[23]] 


TTlh. 


1227 ±4 


1229.5 ± 3.2 


[21] 




3177 


(236 ± 23) 


[24]] 




138.3 ± 1.2 


134.9766 ± 0.0006 


[21] 




1891 ± 15 


1300 ± 100 


[21] 


U 


75.03 ± 0.38 


92.4 ±0.35 


[21] 


rriai 


1114 ±8 


1320 ± 40 


[21] 


fai 


410.9 ±2.4 


433 ± 13 


[21] 




5.54 ±0.02 


6.03 ±0.07 


[21] 


dp'-KTv/ dp'-^TT-K 


-0.128 ±0.002 


-0.13 ±0.02 


[26] 


9pI'-K-kI dp^T^-K 


0.029 ±0.006 


0.028 ±0.02 


[26] 


[r(6i wvr) 


108* 


108 ±9 


[21]] 


\DlS{hx wvr) 


0.15 


0.277 ±0.027 


[21]] 



Table 2. Masses, decay constants, couplings, and decay rates used in our fit. All units (when relevant) 
are MeV. The fit includes all observables listed in the table with the exception of those marked with 
square brackets: data from the bi — )■ wtt decay (F and D/S), and the bi and decay constants. 
The decay constants fb^, /J are evaluated at a scale of 2 GeV. Data from the lattice is shown in 
parentheses. Fit error bars refer to the average degree of accuracy permitted by our numerical analysis, 
and do not include any theory error. 

parameters that naturally appear in the numerical fit of our model are 

along with the remaining independent parameters 55, A, zq, fj,, gigtgb, 929^, and 9z9lgb- Five 
of these, 35, fhq, a, A, and zo, can be compared with analogous parameters in the hard wall 
model: 

"IqH = mquzoH au = ctrz^^ , (4.3) 

together with g^n, Ah, zqh. (The pairs {rhg, a) and {rhgn, an) are analogous in the sense that 
they are the coefficients of the two different powers of the dimensionless parameter ^ = z/zq 
in the function v{z) controlling the tachyon vev, (2.10).) In the hard wall, the parameter 
values (75H = 27r, Ah = 3 were fixed by UV matching, while the global fit of the remaining 
parameters yielded m^n = 2.3 MeV, au = (327 MeV )3, and zqu = (346MeV)-i. These 
correspond to rhqu = 6.65 x 10~^ and cjh = 0.844. 
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0.984 ± 0.008 




±a 


2.44 ±0.03 


rhg/rhqu 


0.74 ±0.04 




9b 


0.31 


a /an 


1.047 ± 0.008 




±9i9b95 


0.381 ± 0.006 


A/Ah 


1.113 ±0.009 




929l 


-3.23 ±0.04 




1.081 ±0.006 




939b95 


77.3 



Table 3. Values of free parameters determined by fit to hadronic data. When relevant, the numbers 
quoted are the ratio between the quantity determined from our fit to that determined in the original 
hard wall model, Model B of reference [1], denoted by the subscript "H". The signs of and gigbga 
are correlated, but the overall sign is not determined. As noted in the text, the quantities gj, and 53 
were not included in the fit, as they only appears in ft,i and fj, and in T and D/S, respectively. The 
value listed in the table is obtained by matching directly to /J and to F. 



It is interesting to note that the parameters of the original hard wall framework tend 
to change very little, even those fixed using perturbative QCD. Meanwhile, the parameter 
/X, related to the scaling dimension of runs significantly from its perturbative value of 

lA'pertl = 1- 

Overall, the results are quite good for the masses and decay constants in the p sector. 
Note that the matching of radial p excitations has improved significantly compared to the 
hard wall predictions. The bi meson is estimated to be far too broad, however. The value of 
939b95 is ^is° unusually large. This indicates that the model has some fundamental failure. 

4.2 Results for the magnetic susceptibility of the quark condensate 

We can use the best-fit parameter values in Table 3 to estimate the value of the magnetic sus- 
ceptibility (3.43). Notice however that the magnetic susceptibility depends on a; using (4.2) 
to convert this to the numerical quantity a leaves us with a dependence on the undetermined 
ratio zq/£. The theoretical prediction also depends on the sign of //, which is not uniquely 
determined by the best fit. We have 

Xu = , , (^)^"' • lim (C^4fe) + C,^(k)) . (4.4) 

The quantity limfc_^o {Cwv{k) + Cvw{k)) can be evaluated numerically for the best-fit pa- 
rameter values and we find its value to be « 0.35. Plugging in the remaining quantities 
yields 

Xu ~ 0.185 X sgn(/x)z2 x (^j « 1.8 GeV-^ x sgn(^) (^-j-j . (4.5) 

The magnetic susceptibility is a quantity of growing relevance for a number of experi- 
mental results. In particular it enters into computations of radiative meson decays [35, 36], 
the muon anomalous magnetic moment [37, 38], and chiral-odd proton distribution functions 
[39, 40]. It has been calculated in QCD using a variety of theoretical tools [30-33] and most 
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recently has been studied using holographic methods [19, 34, 41]. It has also been the subject 
of several lattice QCD computations [27-29]. 

In comparing these results to each other and to our result we find it useful, following [27], 
to write the magnetic susceptibility in terms of a parameter Cu as 

Xu = -c«T7%2 = (2-22 GeV-2) (4.6) 

for Nc = 3. Values of in the literature vary from 0.96 to 1.04 for lattice QCD computations 
[27-29], while ~ 1.4 is a typical value obtained from radiative meson decays [35, 36]. Sum 
rule approaches in QCD which can also be applied in conjunction with ideas of holography 
lead to c„ ~ 2 — 2.15 [32, 34] while the resonance sum rules of [41] require Cu = 4. The only 
other fully holographic computation of c„ uses the model of [7] but with parameter values 
fixed to their UV values as in in [9] and leads to an unacceptably small result, Cu — 0.06 [19]. 

The first thing we learn from comparing our result to the literature is that the diamagnetic 
nature of Xu requires that we choose the sign of to be negative. We can then use the above 
results for Xu to determine a possible range of values for zo/i in our model. We find that 
this ratio varies from roughly 1.8 if we fit to the lattice results to 5 for the radiative decay 
results to 14 and 110 if we try to fit to the QCD sum rule results of [32] or the sum rule 
results of [41] respectively. It is natural to view I/zq as the QCD scale Aqqd and l/£ as a 
renormalization scale Mr in the "conformal window" of this class of AdS/QCD models. A 
ratio Mr/AqcD of 1.8 to 5 is compatible with this point of view while the larger values of 
Mr/AqcD — 14 — 110 seem implausibly large, since they extend the conformal window much 
too far into the ultraviolet. 

Note also that, taking reasonable values for {zq/£) we find a value of Xu in good agreement 
with other determinations, while [19] do not. This can be taken as evidence for our approach, 
in which we fit the model parameters to data rather than demanding that they maintain their 
naive UV values a la [9, 19]. 

5 Conclusion 

We have fully analyzed the extended hard wall framework laid out in [7] with one additional 
interaction term. This extended model incorporates the 6i mesons and a new tower of p states. 
We included all naive dimension 3 QCD operators on the field theory side. We performed a 
new fit for all free parameters in the model, including those of the original framework. The free 
Lagrangian and boundary conditions we used for introducing these new states are essentially 
unique. The interaction terms we included render the model complete up to dimension six 
(according to 5d dimension counting). We succeeded in reproducing much of the observed 
data, including the masses of higher radial p resonances, and an appropriate estimate for the 
magnetic susceptibility. 

However, our rms error is 29.6% zb 2.6%, and our predictions are quite poor for data 
depending most heavily on the bi sector, such as /{,, r(6i — utt) and D/S. The model 
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prefers parameter values such that the Sturm-Liouville weight function in the bi sector varies 
strongly near the IR boundary relative to weight functions in other sectors. This is because 
when §20^ is large and negative as in Table 3, the pole in W(^i,^^ is at z/zq ~ 1.05, close to 
the IR boundary. g2g^ attains a value ten times greater in magnitude than gigbg5, though 

1/2 

9igb95 should dominate by a factor of A'^c • This indicates a significant departure from the 
"conformal" values of (72 = 0, = 1. Though the p sector masses are affected by (72 as 
well, their weight functions depend on it only very mildly, and the original hard wall model 
attained remarkably good fits. 

Our choice of gravity background may underlie this particular failure of the model. As 
discussed in the introduction, the hard wall framework models low-energy QCD with a "con- 
formal window" in which the coupling does not run. This is clearly a crude approximation, 
though one which has proved successful until now. One could check the effect of the back- 
ground quantitatively by comparing our results to the equivalent soft wall model prediction, 
or by introducing more generic (asymptotically AdS) warp factors as polynomials in z, and 
fitting the coefficients to the data. In the end, however, the model reduces to a set of Sturm- 
Liouville problems. The functions appearing in these eigenvalue equations depend on the 
geometry: both directly through the metric, and indirectly through the form of the Xq vev. 
Since the Sturm-Lioville equations conflate these two effects, even our crude analysis indicates 
that the IR geometry is important. 

One should note also that the results we find for the new couplings are not consistent with 
Nc counting. The Nc order of each coupling constant is reproduced in Table 4. Even taking 
into account that Nc = 3, and that each coupling may also contain some order 1 coefficient, 
the hierarchy of the couplings is inconsistent with Nc counting. One should not necessarily 
take this as an indication of a breakdown in the large Nc expansion, but rather a problem 
with 5d dimension counting. For instance, one can increase the 5d dimension, but not the Nc- 
order of any of the interaction terms we have chosen, by sprinkling in appropriately-contracted 
derivatives. We have neglected such terms here. The results indicate that such terms may 
be important, and that 5d dimension counting breaks down as an organizing principle for 
interaction terms. 



95 


6.28 


Nc'" 


9b 


0.31 


Nc'/' 


±9i9b95 


0.38 


Nc'/' 


939b95 


77.3 


Nc'" 


929I 


-3.23 


Nc' 



Table 4. counting and fit values of coupling constants. 

Though we work in a toy model which does not descend from string theory directly, it is 
nevertheless important to note that we have not taken into account all corrections that arise 
at next-to-leading order in a' . We have simply included the first state with mass of order 
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l/^foi' on the vr Regge trajectory. We made this choice because the operator that generates 
this state effectively renders the model complete (a) up to leading order in perturbative QCD 
dimension and (b) up to the mass of the ai. Yet the results of our fit seem to indicate 
that this next-to- leading in 1/\/q7 operator acquires a very large anomalous dimension, so 
the perturbative QCD dimension counting becomes essentially irrelevant. For instance, one 
should also include a' contributions to the background, or the first excitations in l/^/o! on 
the p and oi trajectories. 

In sum, we have produced a model which partially accounts for the finite slope of QCD 
Regge trajectories. The results of our fit (as well as our successful estimate of the magnetic 
susceptibility) suggest that parameters associated with 0{l/^/~Q') states acquire significant 
anomalous dimensions, while those associated with states which appear from the supergravity 
limit of top-down duals do not. One should explore the refinements to the operator content 
and geometry outlined above to confirm these findings. It would also be very interesting to 
understand what effects (either in holography or directly from QCD) specifically suppress 
strong-coupling corrections to these "supergravity- type" operators. 
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A The tachyon vev 

In this appendix we motivate the precise form of (2.10). The complex field Xij is dual to the 
operator gj^(l — 75)<?j = qL,iQR,j, and thus its conjugate, Xij = X*^, is dual to qn^igLj- If 
Xij is the 4d source associated with X, this means that X appears in the QCD generating 
functional as 

lU = ^ exp |i J (fx{ Re {Xij)qiqj + i Im {Xij)qi^^qj + • • • ) | ^ (A.l) 

Thus Re {Xij) is dual to the scalar bilinear Ofj = qiQj while the imaginary part is dual to the 
pseudoscalar. Furthermore we see that Re {X) is naturally identified with the quark mass 
matrix, Re {X) = Mq. 

In QCD there is a spacetime-independent nonzero condensate for the scalar, (ijiqi) = cr. 
We describe this in holography by turning on a background value for Xq = { Re {X)) that is 
independent of the 4d coordinates. The background field should solve the linearized equations 
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of motion, d^^ {^Om^q) + ^{m'^/£'^)Xo = 0, with its value at 
source: 

1 e^-^ 
Me) = ^ • j^M, . 



e determined by the 



(A.2) 



The power of e is determined by the scahng behavior of the non-normahzable mode, and 
the factors of i are inserted on dimensional grounds. This is completely analogous to the 
expressions (3.7). A is related to the tachyon mass mj^ through A = 2 + a /4 + rn^. 



It is convenient to introduce our basis t"" = 1{1,t) and write X = X'^t"', Mq = M^t"". 



The x'^-independent solution with boundary behavior (A.2) is 



X^^iz) 



1 



4-A 



,4-A + CxZ^ 



(A.3) 



Cx is related to the quark vev. Ordinarily it would be fixed by IR boundary conditions and 
these would determine the vev. In the hard wall model we instead treat the vev as input and 
use holography to determine Cx- More precisely we have 



-ij]^Z^r.AM„...] . (A.4) 



Following the philosophy discussed below (3.32), we take Zgrav[Mg] ^ e*'^'^"'!, where S[Mg] 
is the action (2.5) evaluated on the solution (A.3). The leading contribution from (2.5) is 

1 



2 
1 

2£ 



(fx 



(fx 



^X"(e)(5,X") 
^M"M" 

£5-2A 9 Q 



(4 - A)e^-^ + ACxg' 
e4-A(i + c'xe2A-4) 



A-l 



(fx 



((4 - A)e^-2^ + {ACx - (4 - A)Cx) + 0(e'^-')) 



2^6-2A 

which, after subtracting counterterms and taking e — )• 0, leads to^ 



S[Mg 



(fx 



(2A-4^ ■ 

2£6-2A Q Q 



(2A - A)Cx 



-2A 



We can solve this for Cx and plug back into (A.3). As e —J- we find 

1 



ya 



2(A 



(A.5) 



(A.6) 



(A.7) 



Finally, to arrive at (2.10), we specialize to a quark mass matrix that is proportional to 



the identity, Mq 



m, 



,1 = 2mqt^, and a quark vev which is also proportional to the identity. 



a = 2{qt^q) = 2T,^, with the other components vanishing. Then Xq = 2XQt^ is 



X, 







1 

2^ 



2mJ 



A-3^4-A 



+ 



4(A 



1 



(A. 



''in the first step of (A.5) one implicitly assumes that there is no contribution from the IR boundary. 
Equivalently one assumes that there are IR boundary terms in the action which cancel the contribution from 
integration by parts in the bulk; these are the same boundary terms that would play a role in determining Cx ■ 



- 41 - 



B Eliminating auxiliary fields 



In this appendix we outline the steps that take one from the quadratic part of (2.4) to (3.14). 
We decompose the 5d fields into 4d components. Some of the 4d fields will be auxiliary; we 
eliminate them via their equations of motion and plug the result back into the action. 

The decomposition of Am,Vm into 4d fields was discussed around (2.22), (2.23). The 
quadratic actions for Aj^ and {ip,Az) originate entirely from (2.5). They decouple from the 
rest of the fields and are given by 

S^'^[A^] = ^J^J^^"'ftr[d.A^d.A^^-k^Aj^A^ , (B.l) 

S''n^,Az] = y^^iy-^tr[k^\d.^-A.n'-^ . (B.2) 

They encode the tower of J^*^ = l"*"*" /i/oi mesons and = 0~ r//7r's respectively. Our 
extension of the hard wall model does not modify these sectors at quadratic order. The 
equation of motion and boundary conditions for A-^ derived from (B.l) are identical to those 
in [1]. The equations of motion of for the pion sector of [1], described by 5d scalar fields 
(vr, (p), can also be recovered by making the field redefinition dz^p — A^ = dz4>, = (p — ir. This 
field redefinition can be realized as a gauge transformation that sets Az to zero at the price 
of reintroducing pion fluctuations in the phase of X. We prefer to work in a gauge where 
X = Xq is frozen and Az is nontrivial. 

Integration by parts puts (B.l) in the form quoted in the text, with a contribution 
IdM ['"^(a)-^^^^-^^] *o boundary action. For the pion sector we eliminate Az via its 
equation of motion and plug the result back into (B.2), obtaining 

S^'^W] = -j ^£^..I)„tr ||9.^p-A:2M2 + i^M2} , (B.3) 

where 

The unusual momentum dependence in the weight function makes this action awkward to 
work with. It can be shown that solutions are completely regular through the pole in 
but basic properties we expect, such as a complete set of orthogonal wave functions and 
positive-definiteness of the spectrum, are obscured. We can bypass these difficulties by a 
judicious change of field variable. Let 

'^(^' ^) = - 1.2 I ,2 ^zA{k, z) , (B.5) 

and consider the action 

S^'\A\ = I ^£^zt.(.)tr + . (B.6) 
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On shell these two actions are the same: (B.5) maps stationary points of (B.6) to stationary 
points of (B.3). The values at stationary points are also the same. Using the equations of 
motion one can show Jqj^ tr W(^T^)AdzA = — Jgj^ tr [w(^^-^(fdztp]- This is the contribution 
to the boundary action in (3.14). Furthermore (B.6) has a standard form which allows one 
to demonstrate positive definiteness of the spectrum. 

The vector gauge field mixes with the two- form buN at quadratic order due to the gi 
interaction term. The 4d vector b^z is non-dynamical and can be eliminated via its equation 
of motion. To linear order in the fields one finds 

b^,z{k, z)= - Jg^(^)^ ^^i^pg-^^ [\^\bp^{k, z) - U2bpa{-k, z)] + 

- ^^l^^v{z)dzV^ik, z) + 0(cl>2) . (B.7) 

We plug (B.7) back in to yield a quadratic action mixing and the anti-symmetric two-tensor 
bp,v. 

It is convenient to split 6^1^ into its real and imaginary parts, 

bp.v = hp_y + iw^u , (B.8) 

such that h{k,z) = h{—k,z) and similarly for w, and to further split h,w into components 
that are transverse and longitudinal to the momentum k: 

V = ^M*^ + ^lu > 'Wf^iy = wj^u + ■ (B.9) 

The perpendicular and parallel components can be defined in terms of the projectors "P-*-'! 
introduced in (3.34); the only facts we will require are that k^bjj^^^ = and ef_^''"kyb\a = 0. 
Then we find that the remaining piece of the quadratic action splits into two sectors. We 
have the (/i"*", ifH) sector with action 



S^^\h^,w\\] = — I — I dztil sgn{fi)e>^-^P'^CK.dzwl - wldzhj^)+ 



- n2)^"-.l. - ^ (1.1 - - ^) -H^^X^y 

+ ^Af|(,x,^II) , (B.IO) 
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and the (u;-*-, /i", sector with action 

1 r fl'^h /■■^o 



Z ^ \ l/^l ^2 / J 

5af refers to the boundary term of Arutyunov and Frolov in the action (2.7) [7, 13]: 



^2 / 

h JdM 

1 r d^k 



^^9b JdM 



^^9tJ (2vr) 



tr 



(B.12) 



e 



and in (B.IO) and (B.ll) we have restricted it to terms containing the indicated fields only. 

These actions are first order with respect to the h and w fields. On shell they are equiv- 
alent to more standard second order actions that can be obtained by solving the equations of 
motion for the parallel components (say) in terms of the perpendicular ones. To linear order 
in the fields we find 



t^'^ 0(1 ,,\ -U t^'^ ^ per 



e,rdzhj, + 0{<S>^) , (B.13) 



2(1^1 +n2 

which can be substituted back into (B.IO), (B.ll) to get second order actions S^'^'^[h-^] and 

The real transverse anti-symmetric two-tensors , each contain three independent 
components which together comprise the six degrees of freedom in bMN- The transformation 
properties of Imn under charge conjugation and parity identify the normalizable modes of 
/i"*- with the J^*^ = 1"' hi/bi mesons, and the normalizable modes of w-^ with J^*^ = 1 
w/p-mesons. More precisely, the normalizable eigenmodes of the coupled w'^-V system are 
dual to the tower of w/p mesons. The fact that V couples only to w-^ and not is a 
consequence of the fact that the interaction terms (2.16) respect the Z2 x Z2 symmetry of 
charge conjugation and parity. It is convenient to parameterize h-^ and in terms of fields 
7ifj,{k, z), Wfj.{k, z) via 



htAk,z) = -'-^e^rkpH^iKz) , w^,ik,z) = -'-^e^^k^W^ik, z) . (B.15) 
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These fields transform respectively as a 4d axial vector and vector which may be taken 
transverse; the coefficients of their normalizable modes are polarization vectors for the towers 
of and w/yO modes. Plugging (B.13)-(B.15) back into (B.IO), (B.ll) leads to the following 
quadratic actions describing the and w/p sectors: 



dz 



\iA + ^2 



k^2 



1 



\lj\+U2 Z 



- U2) 



+ Saf| 



{h^,w\\) 



(B.16) 



and 



+ 



d^k 

d^k 



(2vr)4 
+ 2sgn(/i) 

- sgn(^) 



dz 
^0 

dz 
z 

d^k 
d'^k k 



trW^ 



z 



k'z 



/^l — U2 J 1^1 — U2 



-M+U2] 



tr {d.V^d.Vf, - fc^v^v^} + 



^0 dz 

Ui 



trV^W,, 



(27r)4 Zo \\n\-U2 ' 



^af|(»x,ii)[W,V] 



(B.17) 



After some integration by parts, the sum of (B.l), (B.6), (B.16), and (B.17) can be put 
in the form (3.14), where the boundary action is 

d'^k 



5'bndry 



(2vr)^ 



tr <^ w^^)Ad,A + W(a)^f,d,A'' + 



1 
24 



1 

ZZq 



sgn(^)A;tii 



sgn{n)kzoUi 



4 2 101/ |2 
2 



+ 



Zo(|/i| - ^2) 



\H\ -U2 



+ 



(B.18) 



A nice feature of this boundary action is that the contribution from z = zq vanishes when 
the IR boundary conditions (3.8) are imposed. Thus, on shell, (B.18) restricts to the z = e 
terms only. 



C Near-boundary analysis 

In this appendix we discuss the asymptotic analysis that leads from the on-shell action at 
finite e, (3.26), (3.30), to the renormalized action in the e — )■ limit, (3.31). Let us begin 
with the pion sector. Examination of the equation of motion, (/c^ — 0^^'^)A = 0, near z = 
shows that the leading behaviors of the non-normalizable and normalizable solution are 

/--=e6-2A(^l + 0(^^min{2,8-2A}^^ ^n.^ ^ ^ ^ (^^ min {2,8-2A}^ _ ^^.1) 
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We have introduced the dimensionless variable ^ = zjzi^ and chosen a convenient normahza- 
tion convention for the solutions: namely the coefficient of the leading power of ^ is equal to 
one. This will be our convention in all sectors. Note that (C.l) does not fully specify the 
non-normalizable solution since we can add a multiple of the normalizable solution without 
changing the asymptotics. This ambiguity would manifest itself in a near boundary analysis 
as a subleading coefficient that is not determined by the differential equation. We define the 
non-normalizable modes so that this coefficient is set to zero. 
The P's and Q's are given by {Pl^' ^P^^^) = (/(°,^ and 



W- _2.8-2A l-^^l^ M> 



Hence 



<3W - '^W^W (6 - 2A)£6-2A ' IJ^O (emin{2,8-2A}) _ q^^Q (e min {2A-4,2}^ 



(C.3) 



We focus on this factor only since the source 5(7^) is 0{e^). Expanding the denominator, we 
see that there can be divergent terms (when A > 3) but that none of these terms involve 
C(^)(/c). In particular, terms coming from the (7(7^) term in the denominator always go as 
a positive power of e; while terms coming from the C(,r) term in the numerator give finite 
contributions. Recalling the expression for the bulk-to-boundary propagator, (3.25), this 
means that only the normalizable component of it contributes in (C.3). 

The set of divergent terms in (C.3) can be determined by a purely asymptotic analysis, 
solving for the coefficients order by order in an asymptotic expansion. It follows from the 
form of the equation of motion that these coefficients will be analytic functions of A:^, and 
therefore the divergences can be eliminated by local counterterms. This results in the pion 
contribution to the renormalized action, (3.31). 

The analysis in the oi and hi sectors is similar. In the ai sector we have 

= l + 0(e'"'°^'''"'^^) , /(a) =C'(l + 0(e"""^'''~'^^)) , (C.4) 

while the Q's and P's are (Q^f', Q^,)) = Uf^', /f,)), and 

PJ^^- = O (^emin{0,6-2A}^ ^ pn.^ = 2 (l + O (^e min {2,8-2A}^ ^ ^ ^^ r^^ 

Again the potentially divergent terms come entirely from {Q'^^')~'^ Pf'^' and can be eliminated 
with local counterterms. There may also be a finite contribution from this piece, but it is 
scheme dependent and may be set to zero as well. The physical finite contribution comes 
from the {Q\;^')~^ P°a) ^^rm. 

For the hi sector we have solutions 

^n.n. ^^-|H^+0('^min{2,8-2A}\\ ^ ^ + ^ (^^ min {2,8-2A}\ \ ^ ((. g^ 
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and Q's and P's, 

=0fe-H+min{2,8-2A}\ ^ = ^ (l + Q min {2,8-2A}\ \ _ (^7) 



ImI 

Q-n. ^^^^0/'^min{2,8-2A}\\ ^ gn^ ^ ^ ^^1+ min {2,8-2A}\ ^ 



^0 

In this case one should keep in mind that the source factors also has e dependence; explicitly, 

^2^2|^t|-2 _ ^ ^ 

^{b)^lS'^^,) = -^^2^^^\'T''""P^,l,,paT'''' , (C.8) 
so we are looking for C'(e^l'^l) terms in the analogous expression to (C.3). The story is the 

)n.n.\ — 1 pn.n 

'(b) ) Ub) 

Hb) ) ""{by 

Finally there is the p sector. Let = (/*, 5*)^ denote the set of four linearly independent 
solutions to (A;2-0(^))i2 = 0. The index z runs over four values, z G {n.n.(zi;), n.n.(?;), n.(it;), n.(f )}. 
The leading small z behavior of these solutions is found to be 

jn.n.;n. ^ ^zfI/^I ^1 + {2,8-2A} 

2sgn(/i)(4 - A)gigbg5mgkz^ 



same. Divergent terms come from (Qf^)') ^P(l)^' and can be removed by counterterms; finite 
terms come from {Qfi))')~^P(h)- 



n.n.;n. 

2sgn(/x)(4-A) gigbQb'mqkzl ^ ^=fImI+6-A /-i I /n / tmin{2A-4,8-2A} 

^(-) |^|(t|/x|+6-A)(^|^|+4-A)^3-a? 



(C.9) 



for the tf-type solutions and 

.5-A 



2 sgn(^) (4 - A)gigbg5mgkzQ (l + O ('^min{2A-4,8-2A}\\ 



((4- A)2 -^2)^3-A ' 

/^mm{2,8-2A}\ ^ (^_^q) 



((6 - A)2 - ^2)^3-A 
^2/i + 0/^min{2,8-2A}\\ ^ 



9{v) 

and 

yn.^ ^ 2sgn(^)(4 - A)ffigbg5mgA:z^~'^ ^g_^ {l + O ('^™ii{2A-4,8-2A} 

for the ?;-type solutions. Note that the relative normalization between /, g for each solution 
is fixed by the equations of motion. 

In order to evaluate the on-shell action, (3.30), we need the P and Q matrices. Recall 
that these have components {P°-^-)ap = P(^) [i?"-'^ ^^)], for example, where a,/3 run over w,v. 
The action of the P's and Q's is given in (3.21), (3.20), and applying them to (C.9)-(C.ll) 
we find 

On.n._ (We)l^l(l + 0(e"^)) 0{e^-^) \ 
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(p) 



(p) 



and 



pn. 

^(P) 



0(eH+4-A) 2(l + 0(e'")) 



(C.13) 
(C.14) 

(C.15) 



where m = min{2,8 — 2A}. 

The matrix sandwiched between the sources in (3.30) can be written as 



M = 



/pn.n.NT 
^■^(P) > 



(C(p))^(P?, 



(P)^ 



cfp)(Q(.))" 



^ ~ (Q{'p)^(p)(Q(p? 



T 



-pn.n. 
^(P) 



[1 



A]-^ \B - Dl 



and the leftmost factor expanded using (1 — A) 
behavior of each of these matrices, finding 



1 + A + 



B 



0(e-l'^l+4-^) 



D 



2(e/zo)l'^IC^, 2a, 



(1 + ©(e'")) + ©(e^*^) 



P(p)^(p)(Qfp)')^ , 

(C.16) 

. We compute the leading 

(C.17) 
(C.18) 



and 








A 


— C-ww 




) + Cwv + C,^ . 0(el^l+6-A) + c,, ■ 0(ei2-2A 


A 


— Cww 


0(e2|Ml+6- 


^) + Cw. ■ 0(elH+i2-2A) ^ . (^(^lH+2) ^ , 


A 


— Cww 


0(£2|Ml+4- 




A 


— Cww 


0(£2|m|+1O 


+ (a. + c^w) • + a. • O(e') . 



(C.19) 

Since we assume A < 4, every term in A goes as some positive power of e. Thus we will only 
need to consider a finite number of terms in the expansion of (1 — A)~^ - enough to soak up 
whatever divergences arise from the rest of the factors in the on-shell action. 

Now consider the e dependence of the source ^(p), which should match the leading asymp- 
totic behavior of i? = (W, V)"^. The leading behavior of the W wavefunction is clearly given 
by the non-normalizable tu-type mode (since we assume A < 4), which behaves as e~l'^L 
The leading asymptotic behavior of the V wavefunction however depends on the quantity 
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6 — A — When this quantity is positive the non-normahzable v-type mode dominates 
with leading behavior e^, but when this quantity is negative, the non-normahzable t(;-type 
mode dominates with behavior e^Ia'I+s-^. if ^e wish to consider general values of |;u| then 
the source in (3.22) should be multiphed by e^Ho-M+a-^) , it follows that finite terms 
in the on-shell action come from terms in M^w that go as e^^'^', terms in M^j^ and My^, that 
go as emin{-|Ml -2|mI+6-A)^ tevuis in My^ that go as emin(o -2|m|+12-2A)^ Terms in the com- 
ponents of M that are dominant to these give divergences, while terms that are subleading 
can be neglected in the limit e — >• 0. 

Consistency of holographic renormalization requires that any divergences be determined 
by purely asymptotic analysis. This will fail if any terms containing components of C are 
divergent. We quickly run into trouble when > 6 — A since M^v for example contains 
a term of the form el^lC^^, which fails to offset the factor of 

^-2M+6-A coming 

from the 

product of sources S'(^„)5(„). When \fi\ < 6 — A there are no such divergences present: B by 
itself gives removable divergences while all terms in A'^B, n > 1, vanish in the e — >• limit. 
Similarly the leading order terms in D contribute finite pieces to the on-shell action while all 
higher order terms in D and A^D vanish. We conclude that the bound 

\n\<6-A (C.20) 

must be imposed. When it is, the p-sector contribution to the quadratic renormalized on-shell 
action is 




(C.21) 



which yields the appropriate terms in (3.31). 



D Mode expansions 

The mode expansions in each sector are as follows. For the 5d fields associated with the pion 
we have 

^{k,z)= • (D.l) 

The 4>l^^ the eigenfunctions of 0^^\ associated with the field A, and the remaining factors 
come from evaluating the relations (3.10) mode by mode. The oi tower is contained in 
via 

Aj^{k,z)=g,Y,^;{k)4:^\z) . (D.2) 
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They, together with the pions, are housed in the 5d axial gauge field, Am, according to (2.22). 
In the bi sector we take 

n 

and finally in the p sector. 

The fields ?^,>V, V are embedded into the 5d fields VM,bMN according to (2.23), (3.11), and 
(3.12). Finally, we recall that we work in a gauge where the remaining 5d field, X, is frozen 
to its vev, X = Xo{z), (2.10). 
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